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DESCRIPTION 


MULTIPLY-COMPLEXED ONE -DIMENSIONAL STRUCTURE, MULTIPLY- 
TWISTED HELIX, MULTIPLY -LOOPED RING STRUCTURE AND 

FUNCTIONAL MATERIAL 


Technical Field 

This invention relates to a multiply- 
complexed one -dimensional structure, multiply- twisted 
helix, multiply- looped ring structure and functional 
material, especially suitable for use as highly 
functional materials based on a novel principle. 
Background Art 

For application of a solid material to 
electronic or optical devices, physical properties of 
the material may restrict its applications. For 
example, in case of using a semiconductor material in a 
liQht emitting device, it will be usable in a device of 
an emission wavelength corresponding to the band gap of 
the material, but some consideration will be necessary 
for changing the emission wavelength. Regarding 
physical properties related to semiconductor bands, 
controls by superlattices have been realized. More 
specifically, by changing the period of a superlattice, 
the bandwidth of its subband can be controlled to 
design an emission wavelength. 

Targeting on controlling many-electron-state 



structures by material designs, the Inventor proposed 
many-body effect engineering by quantum dot -bonded 
structures and has continued theoretical analyses { ( 1 ) 


U.S. 

Patent 

5,430,309; 

(2) 

U.S. Patent 

5,663,571; 

(3) 

u.s. 

Patent 

5,719,407; 

(4) 

U.S. Patent 

5,828,090; 

(5) 

U.S. 

Patent 

5,831,294; 

(6) 

J . Appl . Phys . 76, 



2833(1994); (7) Phys . Rev. B51, 10714(1995); (8) Phys . 
Rev. B51, 11136(1995); (9) J. Appl. Phys. 77, 
5509(1995); (10) Phys. Rev. B53, 6963(1996); (11) Phys. 
Rev. B53, 10141(1996); (12) Appl. Phys. Lett. 68, 
2657(1996); (13) J. Appl. Phys. 80, 3893(1996); (14) J. 
Phys. Soc. Jpn. 65, 3952(1996); (15) Jpn. J. Appl. 

Phys. 36, 638(1997); (16) J. Phys. Soc. Jpn. 66, 
425(1997); (17) J. Appl. Phys. 81, 2693 (1997); (18) 
Physica (Amsterdam) 229B, 146(1997); (19) Physica 
(Amsterdam) 237A, 220(1997); (20) Surf. Sci. 375, 
403(1997); (21) Physica (Amsterdam) 240B, 116(1997); 
(22) Physica (Amsterdam) 240B, 128(1997); (23) Physica 
(Amsterdam) IE, 226(1997); (24) Phys. Rev. Lett. 80, 
572(1998); (25) Jpn. J. Appl. Phys. 37, 863(1998); (26) 
Physica (Amsterdam) 245B, 311(1998); (27) Physica 
(Amsterdam) 235B, 96(1998); (28) Phys. Rev. B59, 
4952(1999); (29) Surf. Sci. 432, 1(1999); (30) 
International Journal of Modern Physics B. Vol. 13, No. 
21, 22, pp. 2689-2703, 1999). For example, realization 
of various correlated electronic systems is expected by 
adjusting a tunneling phenomenon between quantum dots 
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and interaction between electrons in quantum dots. Let 
the tunneling transfer between adjacent quantum dots be 
written as t. Then, if quantum dots are aligned in 
form of a tetragonal lattice, the bandwidth of one 
electron state is Teff=4t. If quantum dots form a one- 
dimensional chain, the bandwidth of one electron state 
is Teff=2t. In case of a three-dimensional quantum dot 
array. Teff=6t. That is, if D is the dimension of a 
quantvim dot array, the bandwidth of one electron state 
has been Teff=2Dt. Here is made a review about half- 
filled (one electron per each quantum dot) Mott 
transition (also called Mott-Hubbard transition or 
Mott metal-insulator transition) . Let the effective 
interaction of electrons within a quantum dot be 
written as Ueff, then the Hubbard gap on the part of the 
Mott insulator is substantially described as A=Ueff- 
Teff, and the Mott transition can be controlled by 
changing Ueff or t. As already proposed, the Mott- 
Hubbard transition can be controlled by adjusting Ueff 
or t, using a field effect, and it is applicable to 
field effect devices (Literatures (5), (6), (11) and 
(14) introduced above). 

On the other hand, reviewing the equation of 
A =Ueff-Teff-2Dt , it will be possible to control Mott- 
Hubbard transition by controlling the dimensionality D 
of the system. For this purpose, the Applicant already 
proposed a fractal-bonded structure that can 
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continuously change the dimensionality, and have 
exhibited that Mott -Hubbard transition is controllable 
by changing the fractal dimensions . 

To enable designing of wider materials, it is 
desired to modify and control the dimension of 
materials by methods different from the fractal theory. 
For example, for the purpose of changing the nature of 
phase transition, it is first conceivable to control 
the number of nearest -neighbor elements among elements 
forming a material. 

On the other hand, here is changed the 
attention to ferromagnetic phase transition taking 
place in the fractal-bonded structure. Of course, 
ferromagnetic materials are one of the most important 
magnetic storage materials. When using z as the number 
of nearest-neighbor atoms, Kb as the Boltzmann 
constant, T as temperature, it is known that 
spontaneous magnetization M in the averaging theory 
describing ferromagnetic phase transition satisfies 

M=Tanh(zM/keT) 

The highest among temperatures T leading to solutions 

of Ms^O of this equation is the critical temperature Tc. 
As readily understood from the equation, Tc is 
proportional to z. When assuming a tetragonal lattice, 
since z=2D, it is expected that the critical 
temperature of ferromagnetic phase transition depends 
on the dimensionality of a material. The Inventor 
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executed more exact Monte Carlo simulation, and showed 
that the critical temperature of ferromagnetic 
transition occurring in a fractal -bonded structure 
could be controlled by the fractal dimensions. 

It is therefore an object of the invention to 
provide a multiply- twisted helix complementary with a 
fractal- shaped material and representing a new physical 
property, and a functional material using the multiply- 
twisted helix. 

A further object of the invention is to 
provide a multiply -looped ring structure complementary 
with a fractal -shaped material and representing a new 
physical property, and a functional material using the 
multiply- looped ring structure. 

A still further object of the invention is to 
provide a multiply- complexed one-dimensional structure 
complementary with a fractal -shaped material and 
representing a new physical property, and a functional 
material using the multiply-complexed one- dimensional 
struc'^re. 

Disclosure of Invention 

The Inventor proposes a multiply-twisted 
helix as one of spatial filler structures. This is 
made by winding a spiral on a spiral structure as a 
base like a chromatin structure that a gene represents, 
and by repeating it to progressively fill a three- 



dimensional space. By adjusting the spiral pitch, the 
spatial filling ratio can be selected, and 
dimensionality of a material, i.e. the number of 
nearest -neighbor elements in this structure can be 
modified. 

In other words , here is proposed a multiply- 
twisted helix in which spirals are made up by using a 
spiral structure as the base and using the spiral 
structure as an element. In this structure including 
hierarchically formed multiple spirals, one -dimensional 
vacancies penetrate the structure to form a structure 
as a porous material. However, by adjusting the turn 
pitch of the spirals, the number of nearest -neighbor 
elements can be changed. According to researches by 
the Inventor, the value of critical inter-electron 
interaction of Mott-Hubbard metal-insulator transition 
in this kind of structure can be controlled by the 
spiral pitch. 

The multiply- twisted helical structure may be 
formed regularly; however, in case a multiply- twisted 
helical structure is actually made, bonding positions 
appearing among spiral layers possibly distribute 
randomly. The degree of the randomness can be new 
freedom of material designs. Taking it into 
consideration, for the purpose of clarifying the effect 
of the random distribution, exact simulation was 
conducted. As a result, introduction of randomness has 



been proved to increase the width of the Mott -Hubbard 
gap and enhance the Mott insulation. Therefore, the 
value of critical inter-electron interaction of Mott- 
Hubbard metal- insulator transition can be controlled 
not only by controlling the degree of randomness of the 
spiral turn pitch but also by controlling the degree of 
randomness regarding inter-layer bonding positions. 

Still in the multiply- twisted helix, there is 
also the inter- layer bonding position as a control 
parameter, in addition to the degree of randomness 
regarding the spiral turn pitch and inter- layer bonding 
positions. That is, by controlling inter-layer bonding 
positions, desired material designs are possible. More 
specifically, freedom of parallel movements of inter- 
layer bonding, that is in other words, simultaneous 
parallel movements of inter-layer bonding positions, 
can be used to control the value of critical inter- 
electron interaction of Mott-Hubbard metal-insulator 
transition occurring in the structure. 

The multiply- twisted helix can be used as a 
magnetic material as well. That is, in the multiply- 
twisted helix, the critical temperature for 
ferromagnetic phase transition to occur can be 
controlled by adjusting the turn pitch. 

The inventor also proposes a multiply- looped 
ring structure as another spatial filler structure, 
which is different in structure from the multiply- 
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twisted helix but similar in effect to same. This can 
be obtained by hierarchically forming rings, using a 
ring as a base. The number of nearest-neighbor 
elements can be changed progressively by adjusting the 
number of elements of low-order rings forming high- 
order hierarchies. Thereby, the spatial filling ratio 
can be established, and the dimensionality of the 
material can be modified. In this multiply -looped ring 
structure, the above discussion is directly applicable 
only if the turn pitch in the multiply- twisted helix is 
replaced by the number of elements. 

The inventor further proposes a multiply- 
complexed one -dimensional structure as a more general 
spatial filler structure that includes both a multiply- 
twisted helix and a multiply-looped ring structure. 

This is obtained by hierarchically forming one- 
dimensional structure systems having a finite curvature 
using a one-dimensional structure system having a 
finite curvature as the base. In this case, the number 
of nearest -neighbor elements can be changed 
progressively by adjusting the curvature of low-order 
one-dimensional structure systems forming high-order 
hierarchies. Thereby, the spatial filling ratio can be 
established, and the dimensionality of the material can 
be modified. In this multiply- complexed one- 
dimensional structure, the above discussion is directly 
applicable only if the turn pitch in the foregoing 
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multiply-twisted helix or the number of elements in the 
foregoing multiply-looped ring structure is replaced by 
the curvature. In this multiply- complexed one- 
dimensional structure, those having a finite twisting 
ratio correspond to multiply-twisted helixes whilst 
those having zero twisting ratios correspond to 
multiply- looped ring structures. 

The present invention has been made as a 
result of progressive studies based on the above 
review . 

That is, to overcome the above- indicated 
problems, according to the first aspect of the 
invention, there is provided a multiply-complexed one- 
dimensional structure having a hierarchical structure 
in which a linear structure as an element of a one- 
dimensional structure having a finite curvature is made 
of a thinner one -dimensional structure having a finite 
curvature, comprising: 

at least two layers of the one- dimensional 
structure^ bonded to each other in at least one site. 

According to the second aspect of the 
invention, there is provided a multiply-complexed one- 
dimensional structure having a hierarchical structure 
in which a linear structure as an element of a one- 
dimensional structure having a finite curvature is made 
of a thinner one-dimensional structure having a finite 
curvature, characterized in: 
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exhibiting a nature regulated by setting a 
curvature in case the one-dimensional structure is made 
of thinner one-dimensional structures. 

According to the third aspect of the 
invention, there is provided a multiply-complexed one- 
dimensional structure having a hierarchical structure 
in which a linear structure as an element of a one- 
dimensional structure having a finite curvature is made 
of a thinner one-dimensional structure having a finite 
curvature, characterized in: 

having a dimensionality regulated by setting 
a curvature in case the one-dimensional structure is 
made of thinner one-dimensional structures. 

According to the fourth aspect of the 
invention, there is provided a multiply-complexed one- 
dimensional structure having a hierarchical structure 
in which a linear structure as an element of a one- 
dimensional structure having a finite curvature is made 
of a thinner one -dimensional structure having a finite 
curvature, having a random potential introduced 
therein, and at least two one-dimensional structures 
bonded in at least one site, characterized in: 

a quantum chaos occurring therein being 
controlled by setting the intensity of the random 
potential, by setting the intensity of layer- to -layer 
bonding, by setting the curvature used when forming the 
one-dimensional structure from thinner one -dimensional 
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structures, or by adding a magnetic impurity. 

Control of the quantum chaos produced 
typically relies on setting a bonding force between 
layers . Addition of magnetic impurities contributes to 
a decrease of the bonding force between layers and good 
control of the quantum chaos. 

According to the fifth aspect of the 
invention, there is provided a functional material 
including in at least a portion thereof a multiply- 
complexed one -dimensional structure having a 
hierarchical structure in which a linear structure as 
an element of a one -dimensional structure having a 
finite curvature is made of thinner one-dimensional 
structures having a finite curvature, characterized in: 
at least two layers of the one- dimensional 
structures being bonded to each other in at least one 
site. 

According to the sixth aspect of the 
invention, there is provided a functional material 
including in at least a portion thereof a multiply- 
complexed one -dimensional structure having a 
hierarchical structure in which a linear structure as 
an element of a one-dimensional structure having a 
finite curvature is made of thinner one-dimensional 
structures haying a finite curvature, characterized in: 
the multiply-complexed one-dimensional 
structure exhibiting a nature regulated by setting the 
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curvature used when the one-dimensional structure is 
made of thinner one-dimensional structures. 

According to the seventh aspect of the 
invention, there is provided a functional material 
including in at least a portion thereof a multiply- 
complexed one -dimensional structure having a 
hierarchical structure in which a linear structure as 
an element of a one -dimensional structure having a 
finite curvature is made of thinner one-dimensional 
structures having a finite curvature, characterized in: 
the multiply- complexed one -dimensional 
structure having a dimensionality regulated by setting 
a curvature in case the one-dimensional structure is 
made of thinner one -dimensional structures. 

According to the eighth aspect of the 
invention, there is provided a multiply- twisted helix 
having a hierarchical structure in which a linear 
structure as an element of a spiral structure is made 
of thinner spiral structures, characterized in: 

at least two layers of spiral structures _ 
being bonded in at least one site. 

According to the ninth aspect of the 
invention, there is provided a multiply- twisted helix 
having a hierarchical structure in which a linear 
structure as an element of a spiral structure is made 
of thinner spiral structures, characterized in: 

exhibiting a nature regulated by setting a 
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turn pitch in case the spiral structure is made of 
thinner spiral structures. 

According to the tenth aspect of the 
invention, there is provided a multiply- twisted helix 
having a hierarchical structure in which a linear 
structure as an element of a spiral structure is made 
of thinner spiral structures, characterized in: 

having a dimensionality regulated by setting 
a turn pitch in case the spiral structure is made of 
thinner spiral structures. 

According to the eleventh aspect of the 
invention, there is provided a multiply- twisted helix 
having a hierarchical structure in which a linear 
structure as an element of a spiral structure is made 
of a thinner spiral structure, having a random 
potential introduced therein, and at least two spiral 
structures bonded in at least one site, characterized 
in: 

a quantum chaos occurring therein being 
controlled by setting the intensity of the random 
potential, by setting the intensity of layer- to- layer 
bonding, by setting the turn pitch used when forming 
the spiral structure from thinner spiral structures , or 
by adding a magnetic impurity. 

Control of the quantum chaos produced 
typically relies on setting a bonding force between 
layers. Addition of magnetic impurities contributes to 
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good control of the quantum chaos. Decreasing the 
bonding force between layers also contributes to good 
control of the quantum chaos. 

According to the twelfth aspect of the 
invention, there is provided a multiply- twisted helix 
having a hierarchical structure in which a linear 
structure as an element of a spiral structure is made 
of a thinner spiral structure, and having at least two 
layers of spiral structures bonded in at least one 
site, characterized in; 

the bonding performance between linear 
structures as elements of the spiral structure being 
controlled by a turn pitch in case of forming the 
spiral structure from thinner spiral structures, by the 
bonding force between the layers, or by a fluctuation 
in the bonding site between at least two layers of 
spiral structures. 

According to the thirteenth aspect of the 
invention, there is provided a functional material 
Including in at least a portion thereof a multiply- 
twisted helix having a hierarchical structure in which 
a linear structure as an element of a spiral structure 
is made of thinner spiral structures, characterized in: 
at least two layers of spiral structures in 
the multiply- twisted helix being bonded in at least one 
site . 

According to the fourteenth aspect of the 
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invention, there is provided a functional material 
including in at least a part thereof a multiply- twisted 
helix having a hierarchical structure in which a linear 
structure as an element of a spiral structure is made 
of thinner spiral structures, characterized in: 

the multiply- twisted helix exhibiting a 
nature regulated by setting a turn pitch produced when 
the spiral structure is made of thinner spiral 
structures . 

According to the fifteenth aspect of the 
invention, there is provided a functional material 
including in at least a part thereof a multiply-twisted 
helix having a hierarchical structure in which a linear 
structure as an element of a spiral structure is made 
of thinner spiral structures, characterized in: 
the multiply- twisted helix having a 
dimensionality regulated by setting a turn pitch in 
case the spiral structure is made of thinner spiral 
structures . 

_ According to the sixteenth aspect of the 

invention, there is provided a multiply-looped ring 
structure having a hierarchical structure in which an 
annular structure as an element of a ring structure is 
made of a thinner ring structure, characterized in: 

at least two layers of ring structures being 
bonded in at least one site. 

According to the seventeenth aspect of the 
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invention, there is provided a multiply-looped ring 
structure having a hierarchical structure in which a 
linear structure as an element of a ring structure is 
made of a thinner ring structure, characterized in: 
exhibiting a nature regulated by setting a 
number of elements in case the ring structure is made 
of thinner ring structures. 

According to the eighteenth aspect of the 
invention, there is provided a multiply- looped ring 
structure having a hierarchical structure in which a 
linear structure as an element of a ring structure is 
made of a thinner ring structure, characterized in: 

having a dimensionality regulated by setting 
a number of elements in case the ring structure is made 
of thinner ring structures . 

According to the nineteenth aspect of the 
invention, there is provided a multiply- looped ring 
structure having a hierarchical structure in which a 
linear structure as an element of a ring structure is 
made of a thinner ring, having a random potential 
introduced therein, and at least two ring structures 
bonded in at least one site, characterized in: 

a quantum chaos occurring therein being 
controlled by setting the intensity of the random 
potential, by setting the intensity of layer- to -layer 
bonding, by setting the number elements used when 
forming the ring structure from thinner ring 
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structures, or by adding a magnetic impurity. 

Control of the quantum chaos produced 
typically relies on setting a bonding force between 
layers. Addition of magnetic impurities contributes to 
a decrease of the bonding force between layers and good 
control of the quantxim chaos. 

According to the twentieth aspect of the 
invention, there is provided a functional material 
including in at least a portion thereof a multiply- 
looped ring structure having a hierarchical structure 
in which a linear structure as an element of a ring 
structure is made of thinner ring structures, 
characterized in: 

at least two layers of the ring structures 
being bonded to each other in at least one site. 

According to the twenty-first aspect of the 
invention, there is provided a functional material 
including in at least a portion thereof a multiply- 
looped ring structure having a hierarchical structure 
in which a linear structure as an element of a ring 
structure is made of thinner ring structures, 
characterized in: 

the multiply- looped ring structure exhibiting 
a nature regulated by setting the number of elements 
used when the ring structure is made of thinner ring 
structures. 


According to the twenty- second aspect of the 
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invention, there is provided a functional material 
including in at least a portion thereof a multiply-loop 
ring structure having a hierarchical structure in which 
a linear structure as an element of a ring structure is 
made of thinner ring structures, characterized in: 

the multiply-looped ring structure having a 
dimensionality regulated by setting a number of 
elements in case the ring structure is made of thinner 
ring structures. 

In the present invention, in the multiply- 
complexed one -dimensional structure, the curvature used 
when a one -dimensional structure in the first layer, 
for example, is made of a thinner one -dimensional 
structure in the second layer lower by one layer than 
the first layer is set to a value different from the 
curvature used when a one-dimensional structure in the 
third layer different from the first layer is made of a 
thinner one-dimensional structure in the fourth layer 
lower by one layer than the third layer. This 
curvature may be set to a different value, depending on 
a difference in position in the one-dimensional 
structure of the same layer. In the multiply- twisted 
helix, the turn pitch used when a spiral structure in 
the first layer, for example, is made of a thinner 
spiral structure in the second layer lower by one layer 
than the first layer is set to a value different from 
the curvature used when a spiral structure in the third 
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layer different from the first layer is made of a 
thinner spiral structure in the fourth layer lower by 
one layer than the third layer. This turn pitch may be 
set to a different value, depending on a difference in 
position in the spiral structure of the same layer. In 
the multiply-looped ring structure, the number of 
elements used when a ring structure in the first layer 
is made of a thinner ring structure in the second layer 
lower than one layer than the first layer is set to a 
value different from the number of elements used when a 
ring structure in the third layer different by one 
layer from the first layer is made of a thinner ring 
structure in the fourth layer lower by one layer than 
the third layer. This number of elements may be set to 
a different value, depending on a difference in 
position in the ring structure of the same layer. 

There may be fluctuation in sites where 
spiral structures, ring structures or one- dimensional 
structures bond between at least two layers. This is 
equivalent to introduction of randomness to sites where 
spiral structures, ring structures or one-dimensional 
structures are bonded between layers. "Fluctuation" 
involves both spatial fluctuation (that can be reworded 
to disturbance or deviation) and temporal fluctuation. 
Any way is employable for introduction of fluctuation. 
For example, it may be introduced to appear in 
predetermined pitches. The fluctuation may be 
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introduced by removing or adding bonds between at least 
two layers of spiral structures, ring structures or 
one -dimensional structures. 

Curvature of the multiply-complexed one- 
dimensional structure, turn pitch of the multiply- 
twisted helix or number of elements of the multiply- 
looped ring structure is made variable under external 
control, for example. In a typical example of the 
multiply-complexed one-dimensional structure, multiply- 
twisted helix or multiply-looped ring structure, the 
one -dimensional structure, spiral structure or ring 
structure is made of a linear formation essentially 
made of atoms or groups of atoms (clusters) as 
elements. To introduce the above-mentioned fluctuation 
into this kind of structure, there is a method of 
inducing random absorption or elimination (surplus 
bonds or lack of bonds) of molecules in the linear 
formation. Since such introduction of a change or 
fluctuation of the curvature, turn pitch of number 
elements can be utilized as a kind of memory function, 
these multiply-complexed one -dimensional structure, 
multiply-twisted helix and multiply- looped ring 
structure can be used as memory devices. 

In the present invention, control of 
curvature in the multiply-complexed one-dimensional 
structure, turn pitch in the multiply- twisted helix or 
number of elements in the multiply -looped ring 
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structure can bring about phase transition, especially 
metal-insulator phase transition or ferromagnetic phase 
transition. Critical value of the phase transition is 
regulated in accordance with the curvature, turn pitch 
or number of elements. This phase transition is 
controlled by control of bonding positions between one- 
dimensional structures, helical structures or ring 
structures of two layers. Specifically, these 
positions are controlled by parallel movements of these 
bonds. Specifically, the multiply-complexed one- 
dimensional structure, multiply- twisted helix or 
multiply-looped ring structure includes, for example, 
metallic phase portions and insulating phase portions . 
The insulating phase portions can change their phase to 
the metallic phase because of their versatility. 
Critical temperature of ferromagnetic phase transition 
can be regulated by the degree of the above-mentioned 
fluctuation. Alternatively, critical temperature for 
ferromagnetic phase transition can be regulated by 
parallel movements of bonds between at least two one- 
dimensional structures, spiral structures or ring 
structures of different layers. 

In a multiply-complexed one -dimensional 
structure, multiply- twisted helix or multiply- looped 
ring structure, including one-dimensional structures, 
helical structures or ring structures of at least two 
different layers, which are bonded, at least, at one 
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position, if this at least one bond itself is made in 
form of a linear structure, then it is possible to 
control various physical phenomena that take place in 
the structure. For example, by setting the force of 
the bond made in form of a linear structure, critical 
temperature for ferromagnetic phase transition can be 
regulated. Alternatively, quantum chaos that may take 
place can be controlled. Further, electron state 
(electron correlation) can be controlled thereby to 
control metal- insulator phase transition. 

In a multiply- complexed one- dimensional 
structure, multiply- twisted helix or multiply- looped 
ring structure, including one-dimensional structures, 
spiral structures or ring structures of at least two 
different layers, which are bonded, at least, at one 
position, if this at least one bond is made via an 
independent element, then it is possible to control 
various physical phenomena that take place in the 
structure. For example, by regulating critical 
-temperature for ferromagnetic phase transition by _ 
making use of the criticality obtained by the 
structure, the structure can exhibit a stable physical 
property against minute structural fluctuations. 
Alternatively, control of the quantum chaos or control 
of metal-insulator phase transition is possible. 
Furthermore, it is possible to control the electron 
state thereby to control metal- insulator phase 
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transition. These structures are not sensitive to 
structural fluctuations and exhibit pandemic physical 
properties, so it is easy to mass-produce a material 
exhibiting a constantly uniform physical property. 

Features of the present aspect of the 
invention common to those of the fourth, 11th and 19th 
aspects of the invention are directly applicable to 
functional materials using any of the multiply- 
complexed one-dimensional structure, multiply-twisted 
helix or multiply- looped ring structure. Features of 
the present aspect of the invention common to those of 
the 12th aspect of the invention are also applicable to 
the multiply- complexed one-dimensional structure when 
the helical structures are replaced by one-dimensional 
structures having a finite curvature, and to the 
multiply- looped ring structure when the helical 
structures are replaced by ring structures. They are 
also applicable to functional materials using the 
multiply- complexed one -dimensional structure, multiply- 
twisted helix or multiply- looped ring structure. 

According to the invention having the above- 
summarized configurations, by appropriately setting 
curvature of the multiply- complexed one-dimensional 
structure, spiral turn pitch of the multiply- twisted 
helix or number of elements in the multiply-looped ring 
structure, it is possible to select a desired value of 
spatial filling ratio, thereby modify or control 
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dimensionality of the material and control physical 
properties that the material exhibit. Especially, for 
example, by controlling the number of nearest -neighbor 
elements in the multiply-complexed one-dimensional 
structure, multiply- twisted helix or multiply-looped 
ring structure, it is possible to modify phase 
transitional natures such as Mott -Hubbard metal- 
insulator phase transition or magnetic phase 
transition. Further, by introducing fluctuations of 
bonding sites between at least two one -dimensional 
structures, spiral structures or ring structures of 
different layers, that is in other words, by 
introducing randomness in bonding sites between one- 
dimensional structures, spiral structures or ring 
structures of different layers and controlling the 
intensity of the randomness, it is possible to adjust 
the insulation performance by inter-electron 
correlation for wider material designs. Furthermore, 
the quantvim chaos can be controlled by appropriately 
setting the potential intensity and the bonding force 
between layers when introducing random potentials , by 
appropriately setting the curvature when one- 
dimensional structures in the multiply-complexed one- 
dimensional structure are made of thin one-dimensional 
structures, by appropriately setting the turn pitch 
when spiral structures in the multiply-twisted helix 
are made of thin helical structures, by appropriately 
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setting the number of elements when ring structure in 
the multiply-looped ring structure are made of thin 
ring structures, or by addition of magnetic impurities. 

Brief Description of Drawings 

Fig. 1 is a rough diagram that shows actual 
bonding in a multiply- twisted helical structure in case 
of N=4; 

Fig. 2 is a rough diagram that schematically 
shows a multiply- twisted helix according to the first 
embodiment of the invention; 

Fig. 3 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=2 in the first embodiment of the invention; 

Fig. 4 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=4 in the first embodiment of the invention; 

Fig. 5 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=6 in the first embodiment of the invention; 

Fig. 6 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=10 in the first embodiment of the invention; 

Fig. 7 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=20 in the first embodiment of the invention; 

Fig. 8 is a rough diagram that shows a 



density of states obtained by numerical calculation in 
case of U=8 in the first embodiment of the invention; 

Fig. 9 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=10 in the first embodiment of the invention; 

Fig. 10 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=12 in the first embodiment of the invention; 

Fig. 11 is a rough diagram that shows a one- 
dimensional structure of a protein; 

Fig. 12 is a rough diagram that shows a 
multiply- twisted spiral as a two-dimensional structure 
using disulfide bonds; 

Fig. 13 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=8 and N=2 in the second embodiment of the 
invention; 

Fig. 14 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=8 and N=4 in the second embodiment of the 
invention; 

Fig. 15 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=8 and N=6 in the second embodiment of the 
invention; 

Fig. 16 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
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case of U=8 and N=10 in the second embodiment of the 
invention; 

Fig. 17 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=8 and N=20 in the second embodiment of the 
invention; 

Fig. 18 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=10 and N=2 in the second embodiment of the 
invention; 

Fig. 19 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=10 and N=4 in the second embodiment of the 
invention; 

Fig. 20 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=10 and N=6 in the second embodiment of the 
invention; 

Fig. 21 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=10 and N=10 in the second embodiment of the 
invention; 

Fig. 22 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=10 and N=20 in the second embodiment of the 
invention; 

Fig. 23 is a rough diagram that shows a 
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density of states obtained by numerical calculation in 
case of U=8 in the second embodiment of the invention; 

Fig. 24 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=10 in the second embodiment of the invention; 

Fig. 25 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=12 in the second embodiment of the invention; 

Fig. 26 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=8 in the third embodiment of the invention; 

Fig. 27 is a rough diagram that shows a 
density of states obtained by nximerical calculation in 
case of U=10 in the third embodiment of the invention; 

Fig. 28 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=12 in the third embodiment of the invention; 

Fig. 29 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation in case of N=2, 3, 4, 6, 10, 20 in the 
fourth embodiment of the invention; 

Fig. 30 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation in case of N=10, 12, 14, 16, 18, 20 in the 
fourth embodiment of the invention; 

Fig. 31 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 



s=l in the fifth embodiment of the invention; 
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Fig. 32 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the fifth embodiment of the invention; 

Fig. 33 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the fifth embodiment of the invention; 

Fig. 34 is a rough diagreim that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the fifth embodiment of the invention; 

Fig. 35 is a rough diagreun that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the fifth embodiment of the invention; 

Fig. 36 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the fifth embodiment of the invention; 

Fig. 37 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the fifth embodiment of the invention; 

_ Fig. 38 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the fifth embodiment of the invention; 

Fig. 39 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the fifth embodiment of the invention; 

Fig. 40 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
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s=l in the fifth embodiment of the invention; 

Fig. 41 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation in case of N=2, 3, 4, 6, 10. 20 in the 
sixth embodiment of the invention; 

Fig. 42 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when R is changed in case of N=4 in the 
sixth embodiment of the invention; 

Fig. 43 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when R is changed in case of N=6 in the 
sixth embodiment of the invention; 

Fig. 44 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when R is changed in case of N=10 in the 
sixth embodiment of the invention; 

Fig. 45 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
Calculation when R is changed in case of N=20 in the 
sixth embodiment of the invention; 

Fig. 46 is a rough diagram that shows changes 
of Q(r) when N is changed in case of R=0 and a=0 in 
the seventh embodiment of the invention; 

Fig. 47 is a rough diagram that shows changes 
of Q(r) when R is changed in case of N=10 and a=0 in 
the seventh embodiment of the invention; 



Fig. 48 is a rough diagram that shows changes 
of Q(r) when a is changed in case of N=10 and R=0 in 
the seventh embodiment of the invention; 

Fig. 49 is a rough diagram that shows a 
coefficient Ci upon optimum approximation of Q(r) by 
Equation 120 in case of R=0 and a=0 in the seventh 
embodiment of the invention; 

Fig. 50 is a rough diagram that shows a 
coefficient C 2 upon optimum approximation of Q(r) by 
Equation 120 in case of R=0 and a=0 in the seventh 
embodiment of the invention; 

Fig. 51 is a rough diagram that shows a 
coefficient Ci upon optimum approximation of Q(r) by 
Equation 120 in case of N=10 and a=0 in the seventh 
embodiment of the invention; 

Fig. 52 is a rough diagram that shows a 
coefficient C 2 upon optimum approximation of Q(r) by 
Equation 120 in case of N=10 and a=0 in the seventh 
embodiment of the invention; 

Fig. 53 is a rough diagram that shows a 
coefficient Ci upon optimum approximation of Q(r) by 
Equation 120 in case of N=10 and R=0 in the seventh 
embodiment of the invention; 

Fig. 54 is a rough diagram that shows a 
coefficient C 2 upon optimum approximation of Q(r) by 
Equation 120 in case of N=10 arid R=0 in the seventh 
embodiment of the invention; 



Fig. 55 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when a is changed in case of N=6 in the 
eighth embodiment of the invention; 

Fig. 56 is a rough diagram that shows 
spontaneous magnetization . obtained by numerical 
calculation when a is changed in case of N=10 in the 
eighth embodiment of the invention; 

Fig. 57 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 

calculation when a is changed in case of N=20 in the 
eighth embodiment of the invention; 

Fig. 58 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the ninth embodiment of the invention; 

Fig. 59 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the ninth embodiment of the invention; 

Fig. 60 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=4 and 
s=l in the ninth embodiment of the invention; 

Fig. 61 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=4 and 
s=l in the ninth embodiment of the invention; 

Fig. 62 is a rough diagram that shows quantum 
level statistics upon changes of N in case of v=2 and 
s=l in the ninth embodiment of the invention; 



Fig. 63 is a rough diagram that shows quantum 
level statistics upon changes of N in case of v=2 and 
s=l in the ninth embodiment of the invention; 

Fig. 64 is a rough diagram that shows quantum 
level statistics upon changes of N in case of v=2 and 
s=0.8 in the ninth embodiment of the invention; 

Fig. 65 is a rough diagram that shows quantum 
level statistics upon changes of N in case of v=2 and 
s=0.8 in the ninth embodiment of the invention; 

Fig. 66 is a rough diagreim that shows quantum 
level statistics upon changes of N in case of v=2 and 
s=0.6 in the ninth embodiment of the invention; 

Fig. 67 is a rough diagram that shows quantum 
level statistics upon changes of N in case of v=2 and 
s=0.6 in the ninth embodiment of the invention; 

Fig. 68 is a rough diagram that schematically 
shows a multiply- looped ring structure according to the 
tenth embodiment of the invention; 

Fig. 69 is a rough diagram that shows a 
density of sJtates obtained by numerical calculation in 
case of U=8 in the tenth embodiment of the invention; 

Fig. 70 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=10 in the tenth embodiment of the invention; 

Fig. 71 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=12 in the tenth embodiment of the invention; 
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Fig. 72 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=6 and 
s=l in the 11th embodiment of the invention; 

Fig. 73 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=6 and 
s=l in the 11th embodiment of the invention; 

Fig. 74 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the 11th embodiment of the invention; 

Fig. 75 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=10 and 
s=l in the 11th embodiment of the invention; 

Fig. 76 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=20 and 
s=l in the 11th embodiment of the invention; 

Fig. 77 is a rough diagram that shows quantum 
level statistics upon changes of v in case of N=20 and 
s=l in the 11th embodiment of the invention; 

Fig. 78 is a rough diagram that shows quantum 
level statistics upon changes of s in case of v=2 and 
L=6 in the 11th embodiment of the invention; 

Fig. 79 is a rough diagram that shows quantum 
level statistics upon changes of s in case of v=2 and 
L=10 in the 11th embodiment of the invention; 

Fig. 80 is a rough diagram that shows quantum 
level statistics upon changes of s in case of v=2 and 
L=20 in the 11th embodiment of the invention; 
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Fig. 81 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=8 in a multiply- looped ring structure of a=0 
as a comparative example in the 12th embodiment of the 
invention; 

Fig. 82 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=10 in a multiply- looped ring structure of a=0 
as a comparative example in the 12th embodiment of the 
invention; 

Fig. 83 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=12 in a multiply- looped ring structure of a=0 
as a comparative example in the 12th embodiment of the 
invention; 

Fig. 84 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=8 in the 12th embodiment of the invention; 

Fig. 85 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=10 in the 12th embodiment of the invention; 

Fig. 86 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of U=12 in the 12th embodiment of the invention; 

Fig. 87 is a rough diagram that shows quantum 
level statistics upon changes of V in case of L=6 and 
s=l in the 13th embodiment of the invention; 
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Fig. 88 is a rough diagram that shows quantum 
level statistics upon changes of V in case of L=6 and 
s=l in the 13th embodiment of the invention; 

Fig. 89 is a rough diagram that shows quantum 
level statistics upon changes of V in case of L=10 and 
s=l in the 13th embodiment of the invention; 

Fig. 90 is a rough diagram that shows quantum 
level statistics upon changes of V in case of L=10 and 
s=l in the 13th embodiment of the invention; 

Fig. 91 is a rough diagram that shows quantum 
level statistics upon changes of V in case of L=20 and 
s=l in the 13th embodiment of the invention; 

Fig. 92 is a rough diagram that shows quantum 
level statistics upon changes of V in case of L=20 and 
s=l in the 13th embodiment of the invention; 

Fig. 93 is a rough diagram that shows quantum 
level statistics upon changes of s in case of L=6 and 
v=2 in the 13th embodiment of the invention; 

Fig. 94 is a rough diagram that shows quantum 
level statistics upon changes of s in case of L=10 and 
v=2 in the 13th embodiment of the invention; 

Fig. 95 is a rough diagram that shows quantum 
level statistics upon changes of s in case of L=20 and 
v=2 in the 13th embodiment of the invention; 

Fig. 96 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when L is changed in case of a=0 in the 
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14th embodiment of the invention; 

Fig. 97 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 

calculation when L is changed in case of a=0.5 in the 
14th embodiment of the invention; 

Fig. 98 is a rough diagram that shows actual 
bonding in a helix-based hierarchy of multiple 
connections according to the 15th embodiment of the 
invention; 

Fig. 99 is a rough diagram that shows changes 
of the average number of nearest -neighbor sites with N 
in a helix-based hierarchy of multiple connections and 
a multiply- twisted helix according to the 15th 
embodiment of the invention; 

Fig. 100 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation in case of N=3. 4, 6. 8, 10. 20 in the 15th 
embodiment of the invention; 

Fig. 101 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=3 in the 
15th embodiment of the invention; 

Fig. 102 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=4 in the 
15th embodiment of the invention; 

Fig. 103 is a rough diagram that shows 



spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=6 in the 
15th embodiment of the invention; 

Fig. 104 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=8 in the 
15th embodiment of the invention; 

Fig. 105 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=10 in the 
15th embodiment of the invention; 

Fig. 106 is a rough diagram that shows 
spontaneous magnetization obtained by n\imerical 
calculation when s is changed in case of N=20 in the 
15th embodiment of the invention; 

Fig. 107 is a rough diagram that shows 
quantum level statistics upon changes of v in case of 
N=10 and s=l in the 16th embodiment of the invention; 

Fig. 108 is a rough diagram that shows 
^quantum level statistics upon changes of v in case of 
N=10 and s=l in the 16th embodiment of the invention; 

Fig. 109 is a rough diagram that shows 
quantum level statistics upon changes of N in case of 
v=2 and s=l in the 16th embodiment of the invention; 

Fig. 110 is a rough diagram that shows 
quantum level statistics upon changes of N in case of 
v=2 and s=l in the 16th embodiment of the invention; 



Fig. Ill is a rough diagram that shows 
quantum level statistics upon changes of s in case of 
N=10 and v=2 in the 16th embodiment of the invention; 

Fig. 112 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=10 and s=l in the 17th embodiment of the 
invention; 

Fig. 113 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of s=l and U=12 in the 17th embodiment of the 
invention ; 

Fig. 114 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=6 and U=6 in the 17th embodiment of the 
invention; 

Fig. 115 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=10 and U=6 in the 17th embodiment of the 
invention; 

Fig. 116 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=20 and U=6 in the 17th embodiment of the 
invention; 

Fig. 117 is a rough diagram that shows actual 
bonding in a critical spiral structure according to the 
18th embodiment of the invention; 

Fig. 118 is a rough diagram that shows 



spontaneous magnetization obtained by numerical 
calculation in case of N=3, 4, 6. 8, 10, 20 in the 18th 
embodiment of the invention; 

Fig. 119 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=3 in the 
18th embodiment of the invention; 

Fig. 120 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=4 in the 
18th embodiment of the invention; 

Fig. 121 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=6 in the 
18th embodiment of the invention; 

Fig. 122 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=8 in the 
18th embodiment of the invention; 

Fig. 123 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=10 in the 
18th embodiment of the invention; 

Fig. 124 is a rough diagram that shows 
spontaneous magnetization obtained by numerical 
calculation when s is changed in case of N=20 in the 
18th embodiment of the invention; 
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Fig. 125 is a rough diagram that shows 
quantum level statistics upon changes of v in case of 
N=10 and s=l in the 19th embodiment of the invention; 

Fig. 126 is a rough diagram that shows 
quantum level statistics upon changes of v in case of 
N=10 and s=l in the 19th embodiment of the invention; 

Fig. 127 is a rough diagram that shows 
quantum level statistics upon changes of N in case of 
v=2 and s=l in the 19th embodiment of the invention; 

Fig. 128 is a rough diagram that shows 
quantum level statistics upon changes of N in case of 
v=2 and s=l in the 19th embodiment of the invention; 

Fig. 129 is a rough diagram that shows 
quantum level statistics upon changes of s in case of 
v=2 and N=10 in the 19th embodiment of the invention; 

Fig. 130 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=10 and s=l in the 20th embodiment of the 
invention; 

Fig. 131 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of s=l and U=6 in the 20th embodiment of the 
invention; 

Fig. 132 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=10 and U=6 in the 20th embodiment of the 
invention; 
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Fig. 133 is a rough diagram that shows a 
density of states obtained by numerical calculation in 
case of N=10 and U=6 in the 20th embodiment of the 
invention; 

Fig. 134 is a rough diagreim that shows a 
density of states obtained by numerical calculation in 
case of N=20 and U=6 in the 20th embodiment of the 
invention; 

Fig. 135 is a perspective view for explaining 
a manufacturing method of a multiply- twisted helix of 
q=2 and N=8 in the 21st embodiment of the invention; 

Fig. 136 is a perspective view for explaining 
a manufacturing method of a multiply-twisted helix of 
q=2 and N=8 in the 21st embodiment of the invention; 

Fig. 137 is a perspective view for explaining 
a manufacturing method of a multiply- twisted helix of 
q=2 and N=8 in the 21st embodiment of the invention; 

Fig. 138 is a perspective view for explaining 
a manufacturing method of a multiply- twisted helix of 
q=2 and N=8 in the 21st embodiment of the invention; 
and 

Fig. 139 is a perspective view for explaining 
a manufacturing method of a multiply- twisted helix of 
q=2 and N=8 in the 21st embodiment of the invention. 

Best Mode for Carrying Out the Invention 

Embodiments of the invention will now be 
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explained below. 
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First explained is a multiply- twisted helix 
according to the first embodiment of the invention. 

An electron system on a multiply- twisted 
spiral in the multiply-twisted helix according to the 
first embodiment is explained below. 

Assuming a one -dimensional lattice, numbers 
are assigned as n=-*-, -1, 0, 1, ". Let the operator for 

generating an electron of a spin a at the p-th lattice 

+■ 

point be Cp„ . Of course, there is the anticommutation 
relation 

( 1 ) 

Here is defined a single -band Hubbard 
Hamiltonian H of the electron system as follows . 


H +i 


Letting electrons be movable only eunong 
neighboring sites, as 

1 when \p-q\ = l 

s when \p-q\ = N 

when vaodi{p,N) = 0andq = p-\-N^ 
when mod{p^N^) = \andq = p + (3) 

s* when mod( p,N^) = 2andq = p + N* 

0 otherwise 

is employed. Assume hereunder that s=l. However, 
mod(a, b) is the remainder as a result of division of a 
by b. Here is made a review about the average number 
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z 


2 2 

2+2 + — + + 
N 


= 2 + 


2N 

N-1 


(4) 

(5) 


Apparently, it can be any value from the value of a 
three-dimensional tetragonal lattice, namely, z=6 when 
N=2, to the value of a two-dimensional tetragonal 
lattice, namely, z=4 when A multiply- twisted 

helical structure is defined by the definition of the 
nearest -neighbor sites- Fig. 1 schematically show how 
the actual bonding is when N=4. Fig. lA is for N pitch. 
Fig. IB is for pitch, and Fig. 1C is for pitch. 
When the structure is folded such that the nearest - 
neighbor sites become spatially closer, the multiply- 
twisted spiral is obtained as shown in Fig. 2. Fig. 2, 
however, embellishes it to provide an easier view. In 
this case, the one-dimensional chain is bonded at two 
right and left sites, and the chain forms spirals of 
the N pitch. Since the spirals form spirals of the 
pitch, the term of enters as a result of transfer of 
electrons between adjacent spirals (Equation 3). Then, 
the spirals fprm spirals of a larger, i.e. N^, pitch. 

Here is defined a pin a electron density 
operator of the j-th site and its sum 

For the purpose of defining a temperature 
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Green's function, here is introduced a grand canonical 
Hamiltonian K = H - fjN where ^ the half 

filled taken here, chemical potential is M=U/2. The 
half-filled grand canonical Hamiltonian can be 
5 expressed as 




iJ,o i 

(6) 


Operators tjj^ , 

/K A 

Jjj,o • ^i,a defined 

beforehand 


as 

t =r^ c r 

(7) 

10 


'i A+ A A+ A 

Jjyi,a “ ^ j,a 

( 8 ) 



u. = c.’.c. .c.\c. , +c. ,c.\ 

* *4 ^4 *4 I,t i4 *4 

( 9 ) 



J A+ A A A + 

( 10 ) 


If the temperature Green function is defined for 

^ A. 

operators A and B given, taking r as imaginary time, 
15 it is as follows. 


( 11 ) 

The on-site Green function 

) = {C;> ;cj,„ } (12) 

is especially important because, when analytic 
20 continuation is conducted as ia>n“*<^+i5 for a small 6 , 

= - (13) 

becomes the local density of states of the site j , and 

D((y) = -^ImGy„(co+t6) (14) 
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becomes the density of states of the system. For later 
numerical calculation of densities of states, 6=0.0001 
will be used. Further let the total number of sites be 
n=10001. 

Imaginary time development of the system is 
obtained by the Heisenberg equation 

■^A{x).[k,A] ( 15 ) 

As the equation of motion of the on-site Green 
function. 



( 16 ) 


is obtained. Then, the approximation shown below is 
introduced, following Gros ((31) C. Gros , Phys. Rev. 
B50, 7295(1994)). If the site p is the nearest- 
neighbor site of the site j, the resolution 

(s.a ^ )(cj,a '^0 ) ( 1? ) 

is introduced as the approximation. This is said to be 
exact in case of infinite-dimensional Bethe lattices, 
but in this case, it is only within approximation. 

Under the approximation, the following equation is 
obtained. 


where 


(io,. 


(18) 



( 19 ) 
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was introduced. To solve the equation, obtained. 


to be analyzed. In case of half -filled 


models, this equation of motion is 

‘■"n {dj,.aCj,a „ - 2/ J {]pJ..oCj,o Ao ] 

Here again, with reference to the Gros logic. 


approximation is introduced. It is the following 
translation. 


( 20 ) 


^ JS. .^ + 

J pJ,-a^j\o 5 ^ j,o 


; c J „ ^ (21) 


Id. c. -c^ ]-*tG id c -c^ 


By executing this translation, the following closed 


equation is obtained. 




(23) 


n j^o 


Here is assumed that there is no dependency on spin. 


That is, assuming 


G, =G., =G , 

1 ;.t jA 


(24) 


the following calculation is executed. 

Densities of States (DOS) 


--y ImG.(o) + i6) 

n ■“ 


that were obtained by numerical calculation are shown 
below. Fig. 3 shows DOS in case of N^2. As shown in 


Fig. 3, when N=2, the density of state D(o)=0) under 
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the Fermi energy to=0 exists in case of U<8, and the 
system behaves as a metal. On the other hand, when 
U=10, a region where DOS has disappeared exists near 
the Fermi energy, and the system behaves as a Mott 
insulator. Fig. 4 shows DOS in case of N=4. As shown 

in Fig. 4, when U=8, D((o=0) becomes substantially 
zero, the system is closer to a Mott insulator. When 
the models of N=6 (Fig. 5). N=10 (Fig. 6) and N=20 
(Fig. 7) are compared, this tendency is enhanced as N 
increases. Regarding the case of changing N while 
fixing U=8, DOS is shown in Fig. 8. From Fig. 8, the 
system is apparently a metal when N=2, and it is easily 
observed that a change to a Mott insulator occurs as N 
increases. Examples of U=10 and U=12 are shown in 
Figs. 9 and 10, respectively. Since a value twice the 
value of to rendering DOS be zero is the gap of a Mott 
insulator (Hubbard gap), it is appreciated that the 
width of the Hubbard gap increases with N. Therefore, 
it has been confirmed that Mott -Hubbard transition can 
be controlled by adjusting N. 

In this manner, when the turn pitch is 
controlled and designed, the system can behave as a 
metal or can behave an insulator. Therefore, a 
material having a plurality of regions having different 
turn pitches can be a superstructure having various 
regions including metallic regions and insulating 
regions, and this enables richer controls and designs 
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of physical properties. For example, if a material is 
designed to have regions of 

metallic/insulating/metallic phases, a device behaving 
as a diode can be realized. Further, if the insulator 
region is made to changeable in phase to a metal under 
external control, the material can behave as a 
transistor. 

A material containing the multiply-twisted 
helix as a part thereof could be made. For example, a 
material in which a structure other than multiply- 
twisted helical structures is added as an electrode to 
a multiply- twisted helical structure will be useful. 

In the foregoing explanation, widths of N, 
and were assumed as turn pitches of multiply- twisted 
spirals; however, this assumption is only for 
simplicity, and there are various other possibilities. 
For example, in a p-th-order spiral, it will be 
possible to select N^+Ap (where Ap is a random number 
of -N^ ^<Ap<N^^). This is realized when a randomness 
regarding positions of bonding points is introduced in 
a high-order spiral. This randomness gives a 
qualitative difference to phase transition of a system. 
It will be explained later in greater detail with 
reference to the second embodiment. 

A method of fabricating such a multiply- 
twisted spiral will, now be explained. If a genetic 
engineering method is used, DNA having a predetermined 
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sequence can be made, a helix is known as one of the 
most popular proteins. This is a spiral structure in 
which four amino acids form one turn. DNA can be made 
such that cysteine, one of one of amino acids and 

having sulfur atoms, is introduced into the a helix. 
When this DNA is copied by mRNA to activate the genes 
to synthesize a protein, its primary structure is as 
shown in Fig. 11. On the other hand, it is known that 
sulfur atoms of cysteine bond one another to form a 
disulfide linkage, and contribute to form high- order 
proteins. It is expected that the multiply- twisted 
spirals will be formed as a secondary structure using 
the disulfide linkage (Fig. 12). 

Next explained is a multiply-twisted helix 
according to the second embodiment of the invention. 

Explained below is an electron system on a 
multiply- twisted spiral in a multiply- twisted helix 
having introduced randomness regarding bonding 
positions between spiral structures of different 
layers . 

Assuming a one -dimensional lattice, numbers 
are assigned as n=-**, -1, 0, I,- *. Let the operator for 

generating an. electron of a spin a at the p-th lattice 

^ + 

point be . Of course, there is the anticommutation 
relation 

' (25) 

Here is defined a single-band Hubbard 
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Hamiltonian H of the electron system as follows. 
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( 26 ) 


Letting electrons be movable only among 
neighboring sites, as 



when |p - = 1 

when |p - ^ I = 

when mod (p,N) = 0 and q = p -i- 

when mod (p,N^) = land q = p + + N^r 

when mod (p,N^) = 2and q = p ■{■ + N^r 


(27) 


0 otherwise 


is employed. Note that mod(a, b) is the remainder as a 
result of division of a by b. r is a random variable 
that satisfies 

-l<r<l (28) 

and the probability that the value of r appears is 
given by 


PRir)^—\r\ (29) 

Average of this distribution is zero, and dispersion is 
the square root of 

r'P^(r)rfr (30) 


1 

1+2//? 


(31) 


When R=l, distribution becomes uniform, and dispersion 
decreases as R decreases. At the limit of R-*0 . 
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randomness disappears. In the following simulation, 
calculation is made for the cases of R=l, R=l/2, R=l/4, 
R=l/8, R=l/16 and R=0. s=l is used below. 

Here is defined a pin o electron density 
5 operator of the j-th site and its sum 

« ( = V « , „ • 

J Zjo J’O 

For the purpose of defining a temperature 
Green's function, here is introduced a grand canonical 
Hamiltonian K = H - fjN where N = • In the half 

10 filled taken here, chemical potential is M=U/2. The 

half-filled grand canonical Hamiltonian can be 
expressed as 

(32) 

/,;,a 


Operators , i2,- and are defined beforehand 

15 as 



i =C^ C +C^ C 

(33) 


^ ^ A+ ^ 

=c c. -c c. 

J j^,o J,q 1,0 lyO ),a 

(34) 

- 


(35) 

20 

J A+ A A A* 

d- = C- c. — c. c. 

i,a i,a /,a i,o i,a 

(36) 


If the temperature Green function is defined for 

operators A and B given, taking t as imaginary time, 
it is as follows. 


{a- b) = -J^ dr(T^A(r)By.‘^"^ ( 37 ) 

25 The on-site Green function 
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<^;.a('wj=(s-.a;cj.c} 

is especially important because, when analytic 
continuation is conducted as icOn-*co+i6 for a small 6, 

^>(") = - + (39) 

becomes the local density of states of the site j , and 

D (a>) = -^ImGj,j(a> + /6) (40) 


becomes the density of states of the system. For later 
nvimerical calculation of densities of states, 6=0.0001 
will be used. Further let the total number of sites be 
n=10001. 

Imaginary time development of the system is 
obtained by the Heisenberg equation 


-~A(t) = [k,a] (41) 

As the equation of motion of the on-site Green 
function, 

is obtained. Then, the approximation shown below is 
introduced, following Gros ((31) C. Gros , Phys . Rev. 
B50 , 7295(1994)). If the site p is the nearest- 
neighbor site of the site j , the resolution 

’^p,a ) ( 43 ) 

is introduced as the approximation. This is said to be 
exact in case of infinite-dimensional Bethe lattices. 
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but in this case, it is only within approximation. 
Under the approximation, the following equation is 
obtained. 




( 44 ) 


where 




( 45 ) 


was introduced. To solve the equation obtained, 
^rfj,.,jC^.o;ct„^has to be analyzed. In case of half -filled 


models, this equation of motion is 






( 46 ) 


Here again, with reference to the Gros logic, 
approximation is introduced. It is the following 


translation. 


(i p,i,-a^i,a » ^j,a ) p, -a j,-a^ j, a > ^],o ) (47) 


; c] o ) ) ( 48 ) 


By executing this translation, the following closed 


equation is obtained. 




( 49 ) 


n ;,a 
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Aspects of DOS under U=8 are shown in Fig. 13 (N=2), 
Fig. 14 (N=4), Fig. 15 (N=6), Fig. 16 (N=10), and Fig. 
17 (N=20). As R increases, randomness also increases. 
Randomness gives almost no influences to metallic 
states (Fig. 13). As N increases and the nature of 
multiply- twisted spirals is introduced, a tail appears 
on DOS under R=0. This tail narrows the Hubbard gap, 
but as the randomness increases, the tail disappears 
and the insulation performance increases, as understood 
from the figure. Aspects of DOS under U=10 are shown 
in Fig. 18 (N=2), Fig. 19 (N=4), Fig. 20 (N=6), Fig. 21 
(N=10), and Fig. 22 (N=20). In Fig. 18, it should be 
remarked that influences of randomness are not so large 
In this case, the system is already a Mott insulator, 
the value of N is small. Therefore, there exists no 
tail structure derived from the multiply-twisted 
helical structure. On the other hand, as N increases, 
a tail is produced, and thereafter disappears due to 
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randomness as understood from the figure. 

It has been confirmed in this manner that, 
when the turn pitch is controlled and designed, the 
system can behave as a metal or can behave an insulator, 
and its insulation performance is enhanced by 
randomness . 

For the purpose of comparison with a system 
in which randomness has been introduced to bonds 
between different layers of the multiply- twisted 
helical structure of different layers, an absolute 
random- bonding system was analyzed. This was directed 
to Mott transition on lattices in which bonds are 
formed absolutely randomly under the constraint that 
each of 10001 sites has at least four nearest -neighbor 
sites. As the average, the number of nearest-neighbor 
sites was determined as 

Z=4 + L/5 (51) 

Therefore, when L=0, it becomes a number of nearest- 
neighbor sites equivalent to that of a two-dimensional 
tetragonal lattice. When L=10, it becomes a number o_f 
nearest -neighbor sites equivalent to that of a three- 
dimensional tetragonal lattice. Aspects of DOS of this 
random- bonding system are shown in Fig. 23 (U=8), Fig. 
24 (U=10) and Fig. 25 (U=12). Changes of DOS are 

smooth, apparently unlike the foregoing multiply- 
twisted helical structure. 

Next explained is an electron system on a 
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multiply- twisted spiral in a multiply- twisted helix 
according to the third embodiment . 

Assuming a one-dimensional lattice, numbers 
are assigned as n=..., -1, 0, 1,.... Let the operator for 
generating an electron of a spin ( at the p-th lattice 
point beCp^ . Of course, there is the anticommutation 
relation 

^p,g^a,p (52) 

Here is defined a single-band Hubbard 

Hamiltonian H of the electron system as follows . 






(53) 


Letting electrons be movable only among 


neighboring sites, as A 


P.9 


ri 


^P.9 ^9.P ' 


S 

2 


when |/7 - = 1 

when |/7 - ^1 = 

when mod(p,N) = 0 and q = p + f{N) 
when mod(/?,7V^^) = 1 and q = p+N^ + f(N^) 
when mod(p,N^) = 2 and q = p + N* + f{N^) 


0 otherwise 


(54) 


is employed. Note that mod(a, b) is the remainder as a 
result of division of a by b. In this definition, an 
arbitrary function f(x) is introduced to provide a 
generalized model. For example, a model given by 

/W=0 (55) 

is the model analyzed in the first embodiment. Whereas, 
a model given by 
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where r is a random variable that satisfies 

-l<r<l (57) 

and having the probability of appearance of r given by 



(58) 


is the model analyzed in the second embodiment. 

A model analyzed in the third embodiment is 
one given by 


f(x) = Qx (59) 

in which Q=-l/2. That is, here is taken a model based 
on the multiply-twisted helix according to the second 
embodiment but having the random variable r fixed to - 
1/2. In other words, it can be considered that inter- 
layer bonds have entirely moved in parallel from the 
system of f(x)=0. 

Here is defined a pin a electron density 
operator of the j-th site 

Aj = V „ . 

J ^ O J>0 

For the purpose of defining a temperature 
Green's function, here is introduced a grand canonical 
Hamiltonian K = H - fxN where ^ . In the half filled 

taken here, chemical potential is n =l3f2. The half- 
filled grand canonical Hamiltonian can be expressed as 




(60) 


i,J,a 


Operators , jj;„ , «, and are defined beforehand as 
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( 61 ) 


^ ^ ^ A + A 

(j4.o =C;.aC/,a+<oCy.a 

h,i,o =<^la^.a -Ko^J,o 
u. =c\c.,c:.c. , +c. ,c\c, ,c. , 

* i,t i,i >.i i,t 

J A+ A A A + 

a =c. c -c. c- 

i,a 1,0 i,a 1,0 1,0 


(62) 

(63) 

(64) 


If the temperature Green function is defined for 
operators A and B given, taking t as imaginary time. 


it is as follows. 


(A;B) = -J^ drlT,A(t)By“^"^ (65) 

The on-site Green function 

Gj,ai^^n) = (cj,a',cla) (66) 

is especially important because, when analytic 
continuation is conducted as ia)n-»0)+i6 for a small 6, 

Dj(co) = - y ImGy „((u + id) (67) 

becomes the local density of states of the site j , and 

D (co) = -^ Im Gy +/6) (68) 

J,o 

becomes the density of states of the system. For later 
numerical calculation of densities of states, 6=0.0001 
will be used. Further let the total number of sites be 
n=10001. 

Imaginary time development of the system is 
obtained by the Heisenberg equation 

-^AC^)-[k,a] (69) 

dr 

As the equation of motion of the on-site Green function. 
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is obtained. Then, the approximation shown below is 
introduced, following Gros ((31) C. Gros, Phys. Rev. 
B50, 7295(1994)). If the site p is the nearest- 

neighbor site of the site j, the resolution 

{^p,a ’ ^j,a ) > ^p,o ^{^i,a » ) ( 71 ) 

is introduced as the approximation. This is said to be 
exact in case of infinite-dimensional Bethe lattices, 
but in this case, it is only within approximation. 
Under the approximation, the following equation is 
obtained. 


(72) 


(73) 


where 

p 

was introduced. To solve the equation obtained, 
(dj_Jj„;c]^„'^has to be analyzed. In case of half-filled 
models, this equation of motion is 




p 

•c^ \ 

■> ho/ 


(74) 

20 

Here again. 

with reference 

to 

the 

Gros logic. 


approximation 

is introduced. 

It 

is 

the following 
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{jpJ,-aCj.aAo)^ ( 75) 

{^j,-a^j,a ’^j,o ) ~*^^p,a J,-o^ j,a ’^],a ^ ( 76 ) 


By executing this translation, the following closed 
equation is obtained. 




[u/lf 


G.. 


io) -t^T. - 2t^T. 

n j,o “ Jt-o 


(77) 


Here is assumed that there is no dependency on spin. 
That is, assuming 


G j = G . . = G . . 

1 ],\ ;.i 


(78) 


the following calculation is executed. In the models 
currently under consideration, there exist parameters 
of the inter- electron interaction U and the turn pitch 
N. For the following calculation, s=t=l is used. 

Let discussion be progressed to a case using 
f(x)=-x/2. Aspects of DOS under N=2, 3, 4, 6, 10, 20 
are shown in Fig. 26 (U=8), Fig. 27 (U=10) and Fig. 28 
(U=12). It will be appreciated that the Mott 
transition is qualitatively controlled by N similarly 
to Fig. 8. However, they are different quantitatively. 
Especially, the Hubbard band tail structure remarkable 
when N is large is small in Fig. 27. Therefore, it is 
possible to modify the electron state in the multiply- 
twisted helix by changing Q in f(x)=Qx. 

Next explained is a multiply-twisted helix 
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according to the fourth embodiment of the invention. 

In the fourth embodiment, ferromagnetic phase 
transition is analyzed by using a Monte Carlo 
simulation to demonstrate that critical temperature for 
5 ferromagnetic phase transition can be controlled by 

adjusting the turn pitch of the multiply- twisted helix. 

Explanation is first made about a spin system 
on a multiply-twisted spiral in the multiply- twisted 
helix. 

10 Assuming a one-dimensional lattice, numbers 

are assigned as n= ", -1, 0, I, "*. Hamiltonian H of 
this spin system is defined as follows. 

( 79 ' 

P.9 

Sp=l or -1 is used as the spin variable, and 
15 an Ising model is handled. As 

1 when |/> - ^1 = 1 
s when \p - q\ = N 

when mod (p,N) = 0 and q = p + 

J p,q = J g,p = • when mod (p,N^) = land q = p + (80) 

s* when mod (p,N^) = 2 and q = p + ~ 

0 otherwise 

is employed. Note that mod(a, b) is the remainder as a 
result of division of a by b. Positional relation 
between spins is determined by J , and the multiply- 
20 twisted helix is defined thereby. Partition function 


62 



( 81 ) 


5 


10 


15 



is introduced, and physical property values are 
calculated at finite temperatures by using statistical 

mechanics. T is the temperature, and p=l/T and kB=l 
are used. Expected value concerning an arbitrary 
function f[Sj) of the spin variable is calculated by 

When n is the total number of sites, expected value of 
spontaneous magnetization is 



(83) 


Hereunder, n=10001 is used. 

Fig. 29 is a result of calculation of 
spontaneous magnetization by using the Metropolis 
method that is one of the Monte Carlo methods. 100000 
is used as the Monte Carlo step. Here are shown values 
obtained by setting N=2, 3, 4, 5, 10, 20 as the turn 
pitch.- The multiply- twisted helix of N=2 exhibits a 
critical temperature substantially equivalent to that 
of a three-dimensional system, but as N increases, the 
critical temperature for spontaneous magnetization to 
disappear gradually decreases. 

The multiply- twisted helixes in which N<10 
exhibit behaviors similar to those of a tetragonal 
lattice in terms of dependency of spontaneous 
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magnetization upon temperature; however, when the 
regions of N>10 are remarked, there are remarkable 
differences in dependency of spontaneous magnetization 
upon temperature. Fig. 30 shows a result of 
calculation of spontaneous magnetization in multiply- 
twisted helixes of N=10, 12, 14, 16, 18, 20 in those 
regions by the same method as used in Fig. 29. As the 
temperature rises, spontaneous magnetization first 
decreases as a differentiable continuous function, and 
thereafter, spontaneous magnetization disappears with a 
differentiable drop. It suggests that this magnetic 
transition is primary phase transition. 

Let a physical interpretation be given. A 
multiply- twisted helix may be considered to be made of 
spirals as one-dimensional elements. Its diameter is 
written as R. Of course, R is a quantity of the same 
degree as the turn pitch N. If the correlation length 
of a spin system is smaller than R, the system is 
considered to be behaving one -dimensionally. Since no 
order parameter is established under one dimension, the 
system is always a disorderly system at finite 
temperatures. However, when the spin correlation 
length exceeds R, spin correlation occurs between 
spiral structures as elements, and the system comes to 
behave as a three-dimensional substance, exhibiting a 
three-dimensional order. That is, it can be 
appreciated that primary phase transition from the one- 
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dimensional disorderly phase to the three-dimensional 
orderly phase . 

Referring to Fig. 30, there are fluctuations 
in temperature for phase transition to occur, and the 
values do not appear to flatly decrease with N; however, 
this is also a feature of primary phase transition. 

Next explained is a multiply- twisted helix 
according to the fifth embodiment of the invention. In 
the fifth embodiment, quantum level statistics will be 
analyzed concerning quantum states in multiply- twisted 
helixes to which random potentials are introduced, 
thereby to demonstrate that quantum chaos can be 
controlled. 

Explained is an electron system on a 
multiply-twisted spiral in a multiply- twisted helix 
according to the fifth embodiment. 

Assuming a one -dimensional lattice, numbers 
are assigned as n= •. -1. 0, 1. •. Let the operator for 
generating a quantum at the p-th lattice point be , 

Of course, let this operator satisfy the 
anticommutation relation 

(84) 

This quantum is a fermion having no free spin freedom. 
This corresponds to analysis of an electron in a solid 
in which spin orbit interaction can be disregarded. 

Here is defined the Hamiltonian H of the 
electron system as follows . 
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( 85 ) 


Letting electrons be movable only among 
neighboring sites, as 

1 when |p - O'! = 1 
s when \p - q^ = N 

when mod (p,N) = 0 and q = p + + f{N) 

when mod (p,N^) = 1 and q == p + + f(N^) (86) 

5 “ when mod (p,N^) = 2 and q = p + + f(N^) 

0 otherwise 

is employed. Note that mod(a, b) is the remainder as a 
result of division of a by b. What is analyzed here is 
a model using 

/(r) = -int^|j ( 87 ) 

where int(x) is the maximum integer that does not 
exceed x. The total number of sites is M, and a 
periodical boundary condition is introduced beforehand 
to p=l, 2, M. The second term of the Hamiltonian 
is the term of the random potential. For each site, 
the random variable 



is generated to form the Hamiltonian. The variable 
breadth v of the random potential is useful as a 
parameter determining the degree of the randomness . 
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Here. N. which determines the turn pitch of multiply- 
twisted spirals, s, which determines the inter- layer 
bonding force, and v, which determines the intensity of 
the random potential, are varied as parameters. 

When eigen energy of the Hamiltonian H is 
, and the eigen vector is |/n) , then 

H\m)^e„\m) (89) 

where m=l, 2, M. 

First, M quantum levels are standardized 
such that their nearest -neighbor level spacing becomes 
1 in average. The density of states of the system is 
defined by 

(90) 

and its staircase function 

X{e) = C ^ drip(r]) (91) 

J — CO 

is calculated. The staircase function obtained is 
converted by using a procedure called unfolding such 
that the density of states becomes constant in average. 
Using the quantum level obtained in this manner, the 
nearest-neighbor level spacing distribution P(s) and 

the Aa statistics of Dyson and Mehta are calculated as 
quantum level statistics. As taught in a literature 
((32) L.E. Reichl, The transition to chaos: in 
conservative classical systems: quantum manifestations 
(Springer, New York, 1992; (33) F. Haake, Quantum 
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Signatures of chaos, ( Springer-Verlag. 1991)), by using 
these statistics, it can be detected whether quantum 
chaos has been generated or not. It is also known that 
a quantum chaotic system is sensitive to perturbation 
5 from outside similarly to the classical chaotic system, 

and analysis of quantum chaos is important as a 
polestar of designs of non-linear materials. 

In case of an integrable system, nearest- 
neighbor level spacing distribution P(s) and A 3 
10 statistics are those of Poisson Distribution 


Pp(s) = e-^ 

(92) 

^ 3 («) = — 

15 

(93) 


In case the system currently reviewed exhibits quantum 
chaos, it becomes GOE (Gaussian orthogonal ensemble) 

15 distribution 


20 


PaoE(s)--f^- 


(94) 


A3(«)=-T 

Jt 


log(2;in)+ y - — — 

8 4 


4-0 


(95) 


where 7 is the Euler's constant. 

Since the multiply- twisted helix analyzed 
here is made up of M=10001 sites, there are M=10001 
eigen states in this quantum system. Among them, the 
following quantum level statistics values were 
calculated on the basis of energy eigen values 
concerning 2001 states from the 4001st to the 6001st 
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states from the ground state. 

Fig. 31 and Fig. 32 show quantum level 
statistics of multiply- twisted helixes of N=10 while 
fixing the inter- layer bonding to s=l and changing the 
ram potential intensity v. Fig. 31 shows P(s) whilst 
Fig. 32 show A 3 statistics values. When v is small, 
the system is in a metallic state, and the quantum 
level statistics is that of a quantum chaotic system of 
GOE. As V increases, the quantum level statistics 
changes toward the Poisson distribution. This is 
typical Anderson transition ((34) Phys. Rev. 109, 1492 
(1958); (35) Phys. Rev. Lett. 42, 783(1979); (36) Rev. 
Mod. Phys. 57. 287(1985)). 

Fig. 33 and Fig. 34 show quantum level 
statistics of multiply- twisted helixes of N=4 while 
fixing the inter-layer bonding to s=l. It is 
recognized that, along with changes from v =2 to v= 22 , 
the system changes from an approximately quantum 
chaotic system to an approximately integrable system. 
What should be noted is the difference the quantum 
level statistics values of v=2 in Figs. 33 and 34 (N=4) 
and the quantum level statistics values of v=2 in Figs. 
31 and 32 (N=10). 

Figs. 35 and 36 shows quantum level 
statistics values obtained when fixing v =2 and s=l and 
changing N. For example, reviewing Aa(n) in detail, it 
is recognized that the quantum level statistics value 
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of v=2 in the multiply-twisted helix of N=4 is larger 
than that in case of N=10. Apparently, as N decreases, 
quantum exhibits a tendency of localization. 

In case of s=l in Figs. 35 and 36, the change 
of the quantum statistics value by the change of N was 
small. This means that the system lacks modulability 
from the viewpoint of the purpose of controlling the 
quantum chaos . Let it tried to change s to improve 
that point. Fig. 37 and Fig. 38 show quantum level 
statistics values obtained when fixing v=2 and s=0.8 
and changing N. Fig. 39 and Fig. 40 show quantum level 
statistics values obtained when fixing v=2 and s=0.6 
and changing N. Similarly to the cases of s=l, quantum 
exhibits the tendency of localization as N decreases; 
however, the breadth of the change is larger than that 
in case of s=l. It has been confirmed that various 
quantum systems can be obtained from substantial 
Poisson distribution to substantial GOE distribution. 

Next explained is a multiply- twisted helix 
according to the sixth embodiment of the invention. In - 
the sixth embodiment , it is demonstrated that 
ferromagnetic transition in a multiply- twisted helix 
can be controlled not only by the number of turns but 
also by randomness of inter- layer bonding. 

Explanation is made about a spin system on a 
multiply- twisted spiral in the multiply- twisted helix 
according to the sixth embodiment. 
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Assuming a one -dimensional lattice, numbers 
are assigned as p=--, -1, 0, I,". Hamiltonian H of 
this spin system is defined as follows. 


H = 




(96) 


Here is employed a three-dimensional vector 
having values at vertices of a regular octahedron. 

That is, Sp=(l,0,0), (-1, 0, 0), (0, 1, 0), (0, -1, 0), 
(0, 0, 1), (0, 0, -1) are allowable as the spin value. 


J =J 

<!.P 


1 when |/j - = 1 

5 when\p-q\ = N 

5^ when mod{p,N) = Q and q = p-\-N^ + f{N) 

■ when mod{p,N^) = \and q = p+N^ + f(N^) 
when mod (j},N^) = 2 and q = p + N* + f(N^) 


(97) 


0 otherwise 


is employed. Note that mod(a, b) is the remainder as a 
result of division of a by b. Positional relation 
between spins is determined by , and the multiply- 
twisted helix is defined thereby. What is analyzed 
here is a model in which randomness is defined in f(x) 
as defined as follows. 

/ (^) = sign {rx) int (|rr|) (98) 

Here is used 


sign (r) = 


I' 1 when r s 0 
I - 1 when r < 0 


(99) 


The term int(x) is the largest integer that does not 
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exceed x. r is the random variable that satisfies 

-l<r<l (100) 

and let the probability that the value of r appears be 
given by 

nw-fkr (101) 

R is the parameter that characterizes the probability 
distribution. Although the average value of the 
distribution is always (r) =0, self multiplication of 
dispersion is given by 

= ( 102 ) 

and dispersion increases as R increases. Under R=0, 
randomness disappears. On the other hand, in case of 
R=l, the distribution becomes most even, and the 
randomness is largest in the model. The random 
variable is generated for each inter-layer bond, and 
randomness is introduced to bonding positions, thereby 
to produce the multiply- twisted helix to be analyzed. 

For the purpose of calculating physical 
values at finite temperatures using statistical 
mechanics , a distribution function 

Z = ^e'^ (103) 

is introduced, where T is the temperature, and P=l/T 
and kB=l are used. Expected value concerning an 
arbitrary function f[s.] of the spin variable is 
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calculated by 


(104) 

When n is the total number of sites, expected value of 
spontaneous magnetization is 

(105) 

Hereunder, n=10001 is used. 

Spontaneous magnetization is calculated by 
using the Metropolis method that is one of the Monte 
Carlo methods. 100000 was used as the Monte Carlo 
steps. First, Fig. 41 shows spontaneous magnetization 
under no randomness of R=0. Here are shown the cases 
setting the turn pitch as N=2, 3, 4, 6, 10, 20 as the 
turn pitch. The multiply- twisted helix of N=2 exhibits 
a critical temperature substantially equivalent to that 
of a three-dimensional system, but as N increases, the 
critical temperature for spontaneous magnetization to 
disappear gradually decreases . This result 
qualitatively coincides with spontaneous magnetization 
in the Ising model in the fourth embodiment in which 
only 1 or -1 is available as the spin. 

Here is made a review about the effect of 
randomness. Spontaneous magnetization in cases using 
R=l/8, 1/4, 1/2 and 1 as the value of R is shown below. 
Fig. 42 shows one under N=4, Fig. 43 shows one under 
N=6, Fig. 44 shows one under N=10, and Fig. 45 shows 
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one under N=20. With each value of N, critical 
temperature changes as R increases. One of its reasons 
might be attenuation of the characterized nature of the 
magnetization curve due to introduction of randomness. 
When N was small, super-period originally unique to the 
multiply- twisted helix did not affect so much, for 
example, in Fig. 42 or Fig. 43, so the magnetization 
curve itself is not different so much from the case of 
a tetragonal lattice. In this case, although the 
critical temperature decreases because of introduction 
of R>0, the shape of the magnetization curve itself 
does not change substantially. On the other hand, when 
N is large, magnetization curves under R=0 in Figs. 44 
and 45, for example, were considerably different from 
that of a tetragonal lattice. In case of R>0, since 
the system progressively changes while losing its 
unique structure, the breadth of the decrease of the 
critical temperature by introduction of randomness is 
large. Through the analysis shown above, it has 
confirmed that critical temperature for ferromagnetic 
transition can be controlled by introduction of 
randomness by R>0. 

Next explained is a multiply- twisted 
structure according to the seventh embodiment of the 
invention. In the seventh embodiment, a simple 
cellular automaton is assumed in the multiply- twisted 
structure, and simulations of its dynamics are 


74 



5 


10 


15 


20 


performed to introduce inter- element spacing. Based on 
the spacing, inter-element bonding property is 
investigated to show that the bonding property can be 
modified. 

Assuming a one- dimensional lattice made up of 
elements, numbers are assigned as p=", -1, 0, I,". A 
multiply- twisted structure is defined by the bonding 
among elements. For the purpose of defining nearest- 
neighbor elements. 



when |/> - = 1 

when \p-q\ = N 

when mod (p,N) = 0and q = p-^-N^ +/(^) 
when mod (jj,N^)= land q = p ■¥ ■¥ f {N'^) 
when mod (p,N^) = 2 and q = p + N* + f(N^) 


(106) 


0 otherwise 


is introduced, and s=l is used hereunder. Assume that 
the p-th element and the q-th element are the nearest- 
neighbor elements with each other when =1. The 
following function is employed as f(x). 

f{x) = sign (rx) int (|rr|) - int (ax) ( 107 ) 


Here was used 


sign(x) =, 


j 1 when r s 0 
I - 1 when r > 0 


(108) 


The term int(x) is the largest integer that does not 
exceed x. r is the random variable that satisfies 

l-<r<l (109) 
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and let the probability that the value of r appears be 
given by 

w-fH"’' (11“) 

R is the parameter that characterizes the probability 
distribution. Although the average value of the 
distribution is always +r,=0, self multiplication of 
dispersion is given by 

(r>)-/<irr=P,(r) = ^ (111) 

and dispersion increases as R increases. Under R=0, 
randomness disappears. On the other hand, in case of 
R=l, the distribution becomes most even, and the 
randomness is largest in the model. The random 
variable is generated for each inter-layer bond, and 
randomness is introduced to bonding positions, thereby 
to produce the multiply- twisted structure to be 
analyzed. a>0 causes parallel movement of the 
distribution. In the multiply- twisted structure, there 
are parameters of the turn pitch N, inter- layer bonds 
distribution width R, inter-layer bonds distribution 
position a. 

Analysis made below is to introduce a certain 
sense of spacing between sites that are not nearest- 
neighbor sites and to measure the inter-site bonding 
property on the basis of the number sites existing at 
the distance from a certain site. For this purpose, a 
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Simple cellular automaton (CA) is introduced. A 
variable of On=l or 0 is introduced to the n-th site. 
Let this variable vary with time t=0, 1, and let it 
written as On(t). 

As the initial condition 


^„( 0 ) = 


1 when n = p 
0 when n ^ p 


( 112 ) 


is used. As the CA dynamics 


10 


15 


is used, where 



0(a:)- 


1 when jc> 0 
0 when jc s 0 


(113) 


(114) 


Through this time development, the time t when the q-th 
site first becomes aq(t)=l is used to define the 
distance Ap, q between the p-th site and the q-th site. 

Let a physical interpretation be given to the 
distance Ap, g. There are a plurality of paths 
connecting the p-th site and the q-th site via the 
nearest -neighbor sites, and the shortest one is the 
above-mentioned distance. In a continuous three- 
dimensional space (x, y, z), it corresponds to a case 
where |x|+|y|+|z| is employed as the distance from the 
origin. In CA mentioned above, propagation to the 
nearest-neighbor site occurs during the time width 1, 
and the time required for Op(0)=l having localized in 
the p-th site to reach the q-th site was measured to be 
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used as the distance. Of course, the distance 


introduced here satisfies the axioms 


P^<1 


(115) 
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A 


P.9 


= A 


9.P 


(116) 
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A„,^A_.+A„. 

P.9 9.»^ 


“ 0 

p>^ 


p = q 


(117) 


By execution of this CA simulation, arbitrary (p, q) 
distance can be determined. 

The quantity remarked in this seventh 
embodiment is the total number (o(r; Pj) of the site q 
at the distance q=r from an arbitrary Pj-th site. 

What is calculated below is the quantity 

Q(r) = ^'^a)(r;pj) (118) 

to be obtained by using a multiply-twisted structure 
made up of M=10001 sites and averaging Ms=1000 samples 
selected as Pj sites from them. It will be readily 
understood that this quantity becomes 

Q(r) = Const. + 4r^ (119) 

in case of a three-dimensional tetragonal lattice. 

Let Q(r) in the multiply- twisted structure 
be calculated. In Fig. 46, Q(r) in case of R=0 and a=0 
was plotted for various N. This system is a multiply- 
twisted structure similar to that taken for analysis of 
Mott transition in the first embodiment. In Fig. 47, 
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Q(r) in case of N=10 and a=0 was plotted for various R 
This system is a multiply- twisted structure similar to 
that taken for analysis of Mott transition in the 
second embodiment. In Fig. 48, Q(r) in case of N=19 
and R=0 was plotted for various a. This system is a 
multiply- twisted structure similar to that taken for 
analysis of Mott transition in the third embodiment. 

In these figures, Q(r) in a region of the level of 
r<10 is described well by quadratic functions of r. 
Hereunder, this is generally developed as 

Q(r) = Cq + C,r^ + Cjr** (120) 

and discussed while using coefficients Ci and C 2 . 

Fig. 49 and Fig. 50 respectively show the 
coefficients Ci and C 2 in the case where Q(r) under 
0<r<10 is optimally approximated by Equation 120 when 
R=0 and a=0 in Fig. 46. Fig. 51 and Fig. 52 
respectively show the coefficients Ci and C 2 in the 
case where Q(r) under 0<r<10 is optimally approximated 
by Equation 120 when N=10 and a=0 in Fig. 47. Fig. 53 
and Fig. 54 respectively show the coefficients Ci and 
C 2 in the case where Q(r) under 0<r<10 is optimally 
approximated by Equation 120 when N=10 and R=0 in Fig. 
48. 

In the cases of Fig. 46 and Fig. 47, the 
Coefficient Ci of O(r^) is modified well by N and R as 
shown in Figs. 49 and 51. As appreciated from Figs. 50 
and 52, since C 2 is very small, Q(r) is described well 
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by quadratic functions. This demonstrates that these 
multiply-twisted structures have three-dimensional 
properties, but they are controlled in number of sites 
in the region separated by the distance r from a 
certain site by modulation of its coefficients . 
Therefore, physical phenomena appearing on the 
structures are controllable. What can be said from 
those results is that, even in the order of r~10, 
evaluation of the average number of the nearest - 
neighbor sites is directly valid and modulability in 
the multiply- twisted helix extends even globally. In 
case of Fig. 48, as apparent from Fig. 53 and Fig. 54, 

there is almost no change in Ci even with changes of a 
unlike the foregoing two examples. On the other hand, 

C 2 varies to a certain degree, and its sign is inverted. 
Therefore, modulability under changes of a of Fig. 48 
results in changes of the coefficient of O(r^). 

Next explained is a multiply- twisted 
structure according to the eighth embodiment of the 
invention. In the eighth embodiment, it is 
demonstrated that ferromagnetic transition in a 
multiply- twisted structure can be controlled by 
parallel movements of inter-layer bonds. 

Assuming a one-dimensional lattice, numbers 
are assigned as p= ", -1, 0, I,-**. Hamiltonian H of 
this spin system is defined as follows. 



( 121 ) 
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H 




p.<? 


Here are employed three-dimensional vectors 
having values at vertices of a regular octahedron. 

That is. Sp=(1.0,0), (-1, 0. 0), (0, 1. 0). (0, -1. 0). 


(0, 0, 1), (0, 0, -1) are allowable as the spin value. 
As J . 




^ P.H *^ 9 ./’ 


1 when |p - ^1 = 1 

s when \p ~(l\- N 

when mod(p,N) = 0 and q = p + + f(N) 

when mod(p,N^) = 1 and q = p + + f(N^) 

s'* when mod(p,N^) = 2 and q = p + 


0 otherwise 

( 122 ) 

is employed. Note that mod(a, b) is the remainder as a 
result of division of a by b. Positional relation 
between spins is determined by J , and the multiply- 
twisted structure is defined thereby. What is analyzed 


here is a model undergoing parallel movements as 
defined by 


/(jc) = -int(ca) (123) 

The term int(x) is the largest integer that does not 
exceed x. In case of a=0, the system exhibits the same 
structure as that of the fourth embodiment. 

For the purpose of calculating physical 
values at finite temperatures using statistical 
mechanics, a distribution function 
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Z = 


r 


is introduced, where T is the temperature, and P=l/T 
and kB=l are used. Expected value concerning an 
arbitrary function f[Sj) of the spin variable is 
calculated by 



(125) 


When n is the total number of sites , expected value of 
spontaneous magnetization is 



(126) 


Hereunder, n=10001 is used. 

Spontaneous magnetization is calculated by 
using the Metropolis method that is one of the Monte 
Carlo methods. Using 100000 as the Monte Carlo steps, 
three different turn pitches, N=6, 10, 20, were 
analyzed. Figs 55 through 57 show results of Monte 
Carlo simulations conducted by using a=0, 0.1, 0.2, 0.3, 
0.4, 0.5. Fig. 55 is the result of N=6, Fig. 56 is the 
result of N=10, and Fig. 57 is the result of N=20. 

Under each of N, critical temperature, i.e^ the 
temperature for spontaneous magnetization to disappear, 
decreases as a increases. This phenomenon would be a 
result of a change of the positional relation between 
nearest -neighbor sites along with an increase of a, 
which makes the spin order fragile. 
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It has been confirmed through the analysis 
that critical temperature for ferromagnetic transition 
can be controlled by parallel movement of inter-layer 
bonds by a>0. 

Next explained is a multiply-twisted helix 
according to the ninth embodiment of the invention. In 
the ninth embodiment, it is demonstrated that quantum 
chaos in a multiply- twisted helix having introduced a 
random potential can be controlled by a random magnetic 
field that can be realized by magnetic impurities. 

A quantum system on a multiply- twisted spiral 
in the multiply- twisted helix according to the ninth 
embodiment will be explained below. 

Assuming a one-dimensional lattice, niambers 
are assigned as p= •, -1, 0, 1, •. Let the operator for 
generating a quantum at the p-th lattice point be . 

Of course, there is the anticommutation relation 

(127) 

This quantum is a fermion having no free spin freedom. 
This corresponds to analysis of an electron in a solid 
in which spin orbit interaction can be disregarded. 

Here is defined the Hamiltonian H of the 
electron system as follows. 

H = t X clc +\ vclc (128) 

P,<l P-<l P 9 ^ P P P ' ' 

P,^ P 

Letting electrons be movable only among 
neighboring sites, as , 
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when |/? - 9 = l| 

when \p-q = N\ 

when mod{p,N) = Q and q=p + N^ + f(N) 

when mod(j?,N^) = l and q = p+N^ + f(N^) (129) 

when mod(p,N^) = 2 and q = p N'* + f(N^) 


0 otherwise 


Note that mod(a, b) is the remainder as a result of 
division of a by b. What is analyzed here is a model 
using 

/(r) = -int^^j (130) 

where int(x) is the maximum integer that does not 
exceed x. tp, q is introduced for the random magnetic 
field. Here is determined as 


where 0p,q=-0q,p is a random real number that satisfies 

O<0p,,<2;r (132) 

Due to hopping to the nearest -neighbor sites, the 
random phase factor 0p,q is re-suffixed, depending the 
sites. If the phase factor is integrated in the loop 
making one turn around a lattice point, a magnetic flux 
passing through the loop is obtained. This means a 
magnetic field is locally introduced to the random 
distribution of O<0p,q<2jt this magnetic field is 
absolutely random .in both intensity and direction, and 
if it is averaged spatially, it becomes a zero magnetic 
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field. 


Using M as the total number of sites, a 
periodical boundary condition is introduced before hand 
to p=l, 2, M. The second term of the Hamiltonian 
is the term of the random potential. For each site, 
the random variable 




(133) 


is generated to form the Hamiltonian. The variable 
breadth v of the random potential is useful as a 
parameter determining the degree of the randomness. 
Here, N, which determines the turn pitch of multiply- 
twisted spirals, s, which determines the inter-layer 
bonding force, and v, which determines the intensity of 
the random potential, are varied as parameters. 

When eigen energy of the Hamiltonian H Is , 
and the eigen vector is , then 

H\m)^e„\m) (134) 

where m=l, 2, M. 

First, M quantum levels are standardized 
such that their nearest -neighbor level spacing becomes 
1 in average. The density of states of the system is 
defined by 


m=l 


and its staircase function 


(135) 
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( 136 ) 


is calculated. The staircase function obtained is 
converted by using a procedure called unfolding such 
that the density of states becomes constant in average. 
Using the quantum level obtained in this manner, the 
nearest -neighbor level spacing distribution P(s) and 

the A 3 statistics of Dyson and Mehta are calculated as 
quantum level statistics. As taught by Literatures 
(32) and (33), by using these statistics, it can be 
detected whether quantum chaos has been generated or 
not. It is also known that a quantum chaotic system is 
sensitive to perturbation from outside similarly to the 
classical chaotic system, and analysis of quantum chaos 
is important as a polestar of designs of non-linear 
materials . 

In case of an integrable system, nearest- 
neighbor level spacing distribution P(s) and A 3 
statistics are those of Poisson Distribution 

Pp(s)=e~^ (137) 

A 3 (n) = ^ (138) 

In case of quantum chaos under a magnetic field, it 
becomes GUE (Gaussian unitary ensemble) distribution 

325^ -4jV 

^ue(^) ~ (139) 

TC 
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= log(2^) + y-7 (140) 

Ut 4 

where T is the Euler's constant. 

Since the multiply- twisted structure analyzed 
here is made up of M=10001 sites, there are M=10001 
5 eigen states in this quantum system. Among them, the 

following quantum level statistics values were 
calculated on the basis of energy eigen values 
concerning 2001 states from the 4001st to the 6001st 
states from the ground state. 

10 Fig. 58 and Fig. 59 show quantum level 

statistics of multiply-twisted structures of N=10 while 
fixing the inter-layer bonding to s=l and changing the 
random potential intensity v. Fig. 58 shows P(s) 
whilst Fig. 59 shows A 3 statistics values. When v is 
15 small, the system is in a metallic state, and the 

quantum level statistics is that of a quantum chaotic 
system of GUE. As v increases, the quantum level 
statistics changes toward the Poisson distribution. 

-This is typical Anderson transition under a magnetic 
20 field ((34) Phys. Rev. 109, 1492 (1958); (35) Phys. Rev. 

Lett. 42, 783(1979); (36) Rev. Mod. Phys. 57. 

287(1985) ) . 

Fig. 60 and Fig. 61 show quantum level 
statistics of multiply-twisted structures of N=4 while 
25 fixing the inter-layer bonding to s=l. It is 

recognized that, along with changes from v =2 to v= 22 . 
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the system changes from an approximately quantum 
chaotic system to an approximately integrable system. 
What should be noted is the difference the quantum 
level statistics values of v=2 in Figs. 60 and 61 (N=4) 
and the quantum level statistics values of v=2 in Figs. 
58 and 59 (N=10) . 

Figs. 62 and 63 shows quantum level 
statistics values obtained when fixing v=2 and s=l and 
changing N. For example, reviewing AaCn) in detail, it 
is recognized that the quantum level statistics value 
of v=2 in the multiply- twisted structure of N=4 is 
larger than that in case of N=10. Apparently, as N 
decreases, quantum exhibits a tendency of localization. 

In case of s=l in Figs. 62 and 63, the change 
of the quantum statistics value by the change of N was 
small. This means that the system lacks modulability 
from the viewpoint of the purpose of controlling the 
quantum chaos . Let it tried to change s to improve 
that point. Fig. 64 and Fig. 65 show quantum level 
statistics values obtained when fixing v=2 and s=0.8 
and changing N. Fig. 66 and Fig. 67 show quantum level 
statistics values obtained when fixing v=2 and s=0.6 
and changing N. Similarly to the cases of s=l, quantum 
exhibits the tendency of localization as N decreases; 
however, the breadth of the change is larger than that 
in case of s=l. As seen from N=2 to N=20, it has been 
confirmed that various quantum systems can be obtained 
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from substantial Poisson distribution to substantial 
GUE distribution. 

Next explained is a multiply- looped ring 
structure according to the tenth embodiment of the 
invention . 


Explained below is an electron system on a 
multiply- looped ring in the multiply- looped ring 
structure according to the tenth embodiment . 

Using a ring having L sites in one cycle, 
multiply-looped rings are defined hierarchically. 
Assuming q hierarchies or layers, the total number of 
sites is L’. The sites are described by a q- 
dimensional cube. 

a =(a„a^,a^-aj (141) 

where ak is an integer that satisfies Isak^sL. 

Let the operator for generating an electron 
of spin a at the a-th lattice point be . Of course, 
there is the anticommutation relation 

Ka4,p|=<5(a-6)5(a-p) (142) 

Here is used 

d(a -fe) = (143) 

j-i 


where 6(a-b) is a delta function defined by 




when a = b 
otherwise 


( 144 ) 


Here is defined a single-band Hubbard 
Hamiltonian H of the electron system as follows . 
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H \X ,c^ c. 

/ j a,b a,o b,o 
a\b~o 


+ U 


y n M , 

fl.t 


Letting electrons be movable only among 
neighboring sites, defines their bonding relation. 
First introduced is the following function. 


^{a,b;k,L) 


fl when -6*1 = 1 or |fl*-6*| = L-l 
lO otherwise 


(146) 


This function has been introduced to provide a 
periodical boundary condition when loops of the period 
L are formed in the k-th layer. 

^(a,b;k,q) = 6(a^ 

(147) 

is a function that takes 1 only when having the same 
suffix, except for the k-th layer. 


q-k 


ri{a,b;k,g)=^Jld{a^^^-j) 




(148) 


has been introduced to prevent bonds from concentrating 
at a particular site. This definition, however, is for 

q>k, and ii(a, b; q, q)=l is affixed beforehand. Using 
these functions 


K,b =K» = ^^^(a,b;k,L)^(a,b;k,q)r](a,b;k,q)s^~’‘ (149) 

is defined. 

Assume hereunder that s=l. However, Let a 
review be made about the average number of the nearest - 
neighbor sites. When the total number of bonds 

=I’+L’+L’-U-. + L^ (150) 


90 



is divided by the total number of sites Lq and then 
doubled, 

2 = 2 + 2^^^ (151) 
is obtained. At the limit of it becomes 

^ ^ L 

z = 2 + 2-^ (152) 

and it is known that it can take any value from the 
value of a three-dimensional tetragonal lattice of z=6 
when L=2 to the value of a two-dimensional tetragonal 
lattice under z=4 when a multiply- looped ring 

structure is defined by using the definition of the 
nearest -neighbor sites. Fig. 68A schematically shows 
the multiply-looped ring structure. Fig. 68B shows a 
ring of L=12 sites. Fig. 68C shows a torus of sites. 
Fig. 68D shows a triple multiply-looped ring structure 
of sites . 

Hereunder, for simplicity, a natural number j 
is defined as 

+ (153) 

In this manner, sites of the multiply- looped ring 
structure can be equated with a one-dimensional lattice 
of IsjsLq. Here is defined a pin a electron density 
operator of the j-th site and its sum 

For the purpose of defining a temperature 
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Green's function, here is introduced a grand canonical 
Hamiltonian K = H-f^ where ^ half 

filled taken here, chemical potential is /x=U/2. The 
half-filled grand canonical Heuniltonian can be 
5 expressed as 


Operators , i2. and are defined beforehand 

as 



A. /v+ A. 

= c c. +C- C- 

jdyO ^y,CT^i,a ''i,a j,a 

(155) 

10 

^ A A 

/ . . = c . C “ c. c . 

J j,i,o '"jyO^'iyO i,a''j,a 

(156) 


u,. =c;,c. ,c.',c. , +c. ,c.',c. ,c.| 

' i.t i.t «,i i.i i,t <,t i.t '.1 

(157) 


d. =c c- — c. c. 

t,a i,a i,a /»a i,<7 

(158) 


If the temperature Green function is defined for 
operators A and B given, taking r as imaginary time. 


15 it is as follows. 

(i; = -J^ (159) 

The on-site Green function 

G,,,(im.) = (c,,„;c’,) ( 160 ) 

is especially important because, when analytic 
20 continuation is conducted as iWn-*w+i5 for a small 6 , 

/);(«) = - + (161) 

becomes the local density of states of the site j, and 

D(cu) = -^ImGy„(co + j6) (162) 
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becomes the density of states of the system. For later 

numerical calculation of densities of states, 5=0.0001 
will be used. 

Imaginary time development of the system is 
obtained by the Heisenberg equation 


^A(z)4k.a] 


(163) 


As the equation of motion of the on-site Green 
function, 

(164) 

is obtained. Then, the approximation shown below is 
introduced, following Gros ((31) C. Gros, Phys. Rev. 
B50, 7295(1994)). If the site p is the nearest- 
neighbor site of the site j, the resolution 

(Cp.c ^ i{Cp,o',clo ;c}_„ ) (165) 

is introduced as the approximation. This is said to be 
exact in case of inf inite- dimensional Bethe lattices, 
but in this case, it is only within approximation. 

Under the approximation, the following equation is 
obtained. 


(ico -t^r b. =1 + — /i. c. -c^ 


(166) 


where 




( 167 ) 


was introduced. To solve the equation obtained. 
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dj_„Cj^;c^j^\ has to be analyzed. In case of half -filled 
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models, this equation of motion is 




■^^^Ki{di,-adp,oAo) 


(168) 


Here again, with reference to the Gros logic, 
approximation is introduced. It is the following 
translation. 


. c ;ct Id- c- ;c^ \ 

\jp,].-a 1 , 0 ’ 1,0 j '^p,-o\ 1,-0 1 , 0 ’ 1 , 0 ! 

Id. c -C^. )^-tG Id. C. ;c^ \ 

\^l,-o'-p,o’^l,o j p,o\^ l,-o'~ 1 , 0 ’'- 1,0 j 


(169) 


(170) 


By executing this translation, the following closed 
equation is obtained. 



= 1 + 


(uHf 


1 ( 0 . 


-rT. -2/T- 

J,o J*-cr 


}’0 


(171) 


Here is assumed that there is no dependency on spin. 
That is, assuming 


G j = G . . = G . . 

J j.t 1,1 


(172) 


the following calculation is executed. 

Values of DOS obtained by numerical 
calculation are shown below. The numerical calculation 
was conducted as to L=3, 4, 5, 6, 10, 20. In case of 
L=3, q=8 is used, and the total number of sites is 
n=6561. Values of the following table were used. 
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respectively. 


L 

q 

n 

3 

8 

6561 

4 

7 

16384 

5 

6 

15625 

6 

5 

7776 

10 

4 

10000 

20 

3 

8000 


Setting U=8 in Fig. 69, U=10 in Fig. 70, U=12 
in fig. 71, and setting L=3, 4, 5, 6, 10, 20 for each 
value of U, numerical calculation was executed. 

In Fig. 69, in case of L=3, the density of 
state D(co=0) under the Fermi energy (0=0 exists, and 
the system behaves as a metal. On the other hand, in 
Figs. 70 and 71, a region where DOS has disappeared 
exists near the Fermi energy, and the system behaves as 
a Mott insulator. When densities of states 
corresponding to different values of L are remarked in 
Fig. 69, in case of L>4, the system has changed to an 
insulator. Therefore, in case of U=8, Mott metal- 
Insulator transition has been confirmed to occur when L 
changes. 

In case of 10, the system behaves as a 

Mott insulator under all analyzed values of L, but 
since the value corresponding to a double of the value 
of (0 rendering DOS be zero is the gap as the Mott 
insulator (Hubbard gap), it is appreciated that the 
width of the Hubbard gap increases with L. Therefore, 
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it has been confirmed that Mott -Hubbard transition can 
be controlled by adjusting L. 

In this manner, when the number of elements L 
of the multiply-looped ring is controlled and designed, 
the system can behave as a metal or can behave an 
insulator. Therefore, a material having a plurality of 
regions having different numbers of elements can be a 
superstructure having various regions including 
metallic regions and insulating regions, and this 
enables richer controls and designs of physical 
properties. For example, if a material is designed to 
have regions of metallic/insulating/metallic phases, a 
device behaving as a diode can be realized. Further, 
if the insulator region is made to changeable in phase 
to a metal under external control, the material can 
behave as a transistor. A material containing the 
multiply-looped ring structure as a part thereof could 
be made. For example, a material in which a structure 
other than the multiply- looped ring structure is added 
as an electrode to a multiply- looped ring structure 
will be useful. 

Next explained is a multiply- looped ring 
structure according to the 11th embodiment of the 
invention. In the 11th embodiment, quantum level 
statistics will be analyzed concerning quantum states 
in the multiply-looped ring structure to which random 
potentials are introduced, thereby to demonstrate that 
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quantum chaos can be controlled. 

Explained below is an electron system on a 
multiply-looped ring in the multiply- looped ring 
structure according to the 11th embodiment. 

Using a ring having L sites in one cycle, 
multiply-looped rings are defined hierarchically. 
Assuming q hierarchies or layers, the total number of 
sites is Lq. The sites are described by a q- 
dimensional cube. 

a = (a^,a^,a^■■■a^) ( 173 ) 

where an is an integer that satisfies l^ak^L. 

Let the operator for generating a quantum at 
the a-th lattice point be Cj . Of course, let the 
operator satisfy the anticommutation relation 

-^) ( 174 ) 

This quantum is a fermion having no free spin freedom. 
This corresponds to analysis of an electron in a solid 
in which spin orbit interaction can be disregarded. 
Here is used 




where 6(a-b) 


is a delta function defined by 

1 when a ==b 
0 otherwise 




( 175 ) 


( 176 ) 


Letting electrons be movable only among 
neighboring sites , determines their bonding 
relations. First introduced is the following function. 
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1 when b^--a^-±\ 
\a,b;k,L) = . 1 when b^-a^ = ±(1 - L) 
0 otherwise 


( 177 ) 


This function has been introduced to provide a 
periodical boundary condition and bonds between 
adjacent loops in the k-th layer. Next introduced is 




(178) 


This is a function that takes 1 only when having the 
same suffix concerning layers beyond the k-th layer. 
Furthermore, by using 

Q{a,b;k,q)=\ ^ Jr\k., ji 'i {119) 

I 1 when q^k-\-\ 


R^'\x,a,L) = 5 mod(r-l± — ,L) 
\ 2 


also introduced is 


Using these functions. 


(180) 


^^%Mk,q) = R^%,,,;a,L)R^^-\b,^,;a,L)Q(aMk,^ ( 181 ) 


a«=+,- jT=1 

is defined. When placing a=0, the above definition 
results in the same as that already discussed in the 
tenth embodiment. When a>0, the same effects as 
parallel movements of inter- layer bonds in the 
multiply- twisted helix is expected. Regarding its 
effects to Mott metal-insulator transition, detailed 
explanation will be made later in the second embodiment, 
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Hereunder, using a=0.5, calculation will be conducted. 

Generally, a natural number j is defined as 

7 (183) 

In this manner, sites of the multiply-looped ring 
5 structure can be equated with a one- dimensional lattice 

of l:sjsLq. The Hamiltonian H of this quantum system is 
defined as 

H 2 ( 184 ) 

a,& a 

The second term of the Hamiltonian is the term of the 
10 random potential. For each site, the random variable 

--<v,<- (185) 

is generated to form the Hamiltonian. The variable 
breadth v of the random potential is useful as a 
parameter determining the degree of the randomness. 

15 When eigen energy of the Hamiltonian H ±s , 

and the eigen vector is \m) , then 

H\m)^eJ^m) (186) 

where m=l, 2, ••, n. 

First, n quantum levels are standardized 
20 such that their nearest -neighbor level spacing becomes 

1 in average. The density of states of the system is 
defined by 

( 187 ) 
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and its staircase function 


^(e)=fjrip{r]) (188) 

is calculated. The staircase function obtained is 
converted by using a procedure called unfolding such 
that the density of states becomes constant in average. 
Using the quantum level obtained in this manner, the 
nearest -neighbor level spacing distribution P(s) and 
the A 3 statistics of Dyson and Mehta are calculated as 
quantum level statistics . As taught by Literatures 
(32) and (33), by using these statistics, it can be 
detected whether quantum chaos has been generated or 
not. It is also known that a quantum chaotic system is 
sensitive to perturbation from outside similarly to the 
classical chaotic system, and analysis of quantum chaos 
is important as a poles tar of designs of non-linear 
materials . 

In case of an integrable system, nearest - 
neighbor level spacing distribution P(s) and A 3 
statistics are those of Poisson Distribution 

Pp{s)^e-^ (189) 

A3(n) = ^ (190) 

In case the system discussed here exhibits quantum 
chaos , it becomes GOE distribution 

^G0£(^) = Y^ (191) 
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(192) 


^ 3 («) = -T 
7T 


\og(2m) + y - 


n 



where 7 is the Euler's constant. 

In this section, L, which determines the turn 
pitch of the multiply- looped ring structure, s, which 
determines the inter-layer bonding force, and v, which 
determines the intensity of the random potential, are 
varied as parameters. When L=6, 10, 20 is used as the 
turn pitch, the values 


L 

q 

n 

6 

5 

7776 

10 

4 

10000 

20 

3 

8000 


were used as the total number of sites. Statistics 
values are calculated by using 2001 quantum levels, in 
total, taking 1000 states on each side from the band 
center among n states used for each L. 

Fig. 72 and Fig. 73 show quantum level 
statistics of multiply- looped ring structures of L=6 
while fixing the inter-layer bonding to s=l and 
changing the random potential intensity v. Fig. 74 and 
Fig. 75 show quantum level statistics of multiply- 
looped ring structures of L=10 while fixing the inter- 
layer bonding to s=l and changing the random potential 
intensity v. Fig. 76 and Fig, 77 show quantum level 
statistics of multiply-looped ring structures of L=20 
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while fixing the inter-layer bonding to s=l and 
changing the random potential intensity v. Figs. 72, 
74 and 76 show P(s) whilst Fig. 73, 75 and 77 show A 3 
statistics values. In each of Figs. 72 through 77, 
when V is small, the system is in a metallic state, and 
the quantum level statistics is that of a quantum 
chaotic system of GOE. As v increases, the quantum 
level statistics changes toward the Poisson 
distribution. This is typical Anderson transition 
((34) Phys. Rev. 109, 1492 (1958); (35) Phys . Rev. Lett. 
42, 783(1979); (36) Rev. Mod. Phys. 57. 287(1985)). 

Through a careful review of the case of v=2, 
it is confirmed that the chaotic tendency is enhanced 
as L increases. On the other hand, in case of v=22, as 
L increases , the system asymptotically changes toward 
the Poisson distribution. That is, as L increases, the 
level statistics values of this quantum system become 
more sensitive to changes of v. 

Similarly to the multiply- twisted helix, here 
is explained is a result of analysis of a system in 
which the inter- layer bonding intensity s is decreased. 

While fixing v=2. Fig. 78 shows A 3 statistics values of 
multiply- looped ring structures of L= 6 , Fig. 79 shows 

A 3 statistics values of multiply- looped ring structures 
of L=10, and Fig. 80 shows A 3 statistics values of 
multiply-looped ring structures of L=20 . When s 
decreases from 1 , the inter- layer bonding becomes 
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weaker, and the spatial dimension reduces. Therefore, 
it results in enhancing localization. It is 
appreciated from Figs. 78 through 80 that the system 
changes from a quantum chaotic sate to a considerably 
localized state along with changes from s=l to s=0.4. 
The degree of the changes becomes larger as L becomes 
smaller. This is because, when L is small, because of 
a small diameter of the basic double ring, this is 
substantially a one-dimensional tube, but in contrast, 
when L is large, because of a large diameter of the 
basic double ring, this itself forms a two-dimensional 
sheet . 

Next explained is a multiply- looped ring 
structure according to the 12th embodiment of the 
invention. In the 12th embodiment, parallel movement 
of inter-layer bonds in multiply-looped ring structures 
are introduced. Thereby, it is demonstrated that, by 
the use of this freedom, the nature of Mott metal- 
insulator transition occurring in multiply- looped ring 
structures can be controlled and inter-electron 
correlation intensity can be designed. 

Explained below is an electron system on a 
multiply- looped ring in the multiply- looped ring 
structure according to the 12th embodiment. 

Using a ring having L sites in one cycle, 
multiply- looped rings are defined hierarchically. 
Assuming q hierarchies or layers, the total number of 
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sites is Lq. The sites are described by a q- 
dimensional cube. 


a^{a„a^,aya^) 


(193) 


Where an is an integer that satisfies IrsaK^L. 

Let the operator for generating an electron 
of spin a at the a-th lattice point be . Of course, 
there is the anticommutation relation 

p) (194) 

Here is used 




(195) 


where 6(a-b) is a delta function defined by 


fl when a=b 
otherwise 


(196) 


Here is defined a single-band Hubbard 
Hamiltonian H of the electron system as follows . 

H =t \ ^cl _ +U \ h .h . (197) 


a,b,a a 


Letting electrons be movable only among 
neighboring sites, defines their bonding relation. 
First introduced is the following function. 




1 when ~ -1 

. 1 when b^-a^ - -(!“!') 
0 otherwise 


(198) 


This function has been introduced to provide a 
periodical boundary condition when loops of the period 
L are formed in the k-th layer. Next, 
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(199) 


^iaMk)-U6(aj-bj) 

;=1 

is a function that takes 1 only when having the same 
suffix concerning layers beyond the k-th layer. 
Furthermore, by using 


1 when q^k-^l 


( 200 ) 


K^\x-,a,L) = (5 1 mod(r - 1 ± — ,L) 


( 201 ) 


also introduced is 


= R^^\a,,,;a,L)R^^\b,^,;a,L)Q(a,b-Xq) ( 202 ) 


Using these functions. 


= S \E‘’(aMk,L)^^%;k,qma,b;k)s‘'-’‘ ( 203 ) 


a=+,- 


is defined. Hereunder, s=l is assumed. 

For simplicity, in the subsequent procedures, 
a natural number j is defined as 


1 




(204) 


Jn this manner, sites of the multiply-looped ring 
structure can be equated with a one -dimensional lattice 


of ls 3 <Lq. Now defined are the operator of density of 


electrons with spin o in the j-th site 


and the sum thereof 




", - “ 


(205) 


(206) 


For the purpose of defining a temperature 
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Green's function, here is introduced a grand canonical 
Hamiltonian 



K = H-fjN 

where 

(207) 

5 

II 

M 

js> 

(208) 


In the half filled taken here, chemical 

potential is U 


=U/2. The half -filled grand canonical 
be expressed as 

Hamiltonian can 


iTjTcf 1 

0 (209) 

10 

Operators tjj^ , , u. and are defined beforehand 

as 


^ A A+ A 

^ j,i,cr "" ^j\o^iyO “ j,o 

A A+ A A f A A A+ A A + 

(210) 


(211) 


(212) 

15 

d- — c — c 

(213) 


If the temperature Green function is defined for 


operators A and B given, taking T as imaginary time, 
it is as follows. 

(214) 

20 The on-site Green function 

) = (S.0 ; ) (215) 

is important for calculation of densities of states. 

Imaginary time development of the system is 
obtained by the Heisenberg equation 
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^^(t) . [k,a] 

dr 


(216) 


As the equation of motion of the on-site Green function, 

= (217) 

p.i ^ 

is obtained. Then, the approximation shown below is 
introduced, following Gros ((31) C. Gros, Phys. Rev. 
B50, 7295(1994)). If the site p is the nearest- 


neighbor site of the site j , the resolution 




(218) 


is introduced as the approximation. This is said to be 
exact in case of infinite-dimensional Bethe lattices, 
but in this case, it is only within approximation. 
Under the approximation, the following equation is 


obtained. 


(to. + (219) 


where 


r,:. - 


( 220 ) 


was introduced. To solve the equation obtained. 


dj_aCj„’,Cj^) has to be analyzed. In case of half-filled 


models, this equation of motion is 


;cl. ) - |c,„ - 2,^ 
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( 221 ) 


Here again, with reference to the Gros logic, 
approximation is introduced. It is the following 


translation. 


/; c ;ct Id- c ;c^ \ (222) 

\J Pyjy~o j,o j p,-o\"j,-a j,a^ j,ol 


Id- c ;ct \—^tG Id- c- ;ct 

j,-a p^o^ jyO / p,cr\ ;,-a ;,a’ ;,a 


(223) 


By executing this translation, the following closed 


equation is obtained. 




(224) 


n JyO 


Here is assumed that there is no dependency on spin. 


That is, assuming 


G.{ico„) = G .^ (iwj = G. , (ia>„) 


(225) 


the following calculation is executed. When analytic 
continuation is conducted as icOn-^co+i <5 for a small 6 , 


Dj{(o) = - Im Gy((W + id) 


(226) 


becomes the local density of states of the site j, and 


n pi 


(227) 


becomes the density of states of spin a electrons . 

Values of DOS obtained by numerical 
calculation while fixing t=l are shown below. 
Regarding L=3, 4, 6, 8, 10, 20, values in the following 
table are used as the number of layers q and the total 
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number of sites n=L'*. 


L 

q 

n 

3 

8 

6561 

4 

6 

4096 

6 

5 

7776 

8 

4 

4096 

10 

4 

10000 

20 

3 

8000 


First, the case of a=0 (tenth embodiment) is 
reviewed again. Numerical calculation was executed 
setting U=8 in Fig. 81, U=10 in Fig. 82, U=12 in Fig. 
83, and selecting L=3, 4, 6, 8. 10, 20 for each value 
of U. 

In Fig. 81, in case of L=3, the density of 
state D(to=0) under the Fermi energy oi=0 exists, and 
the electron system behaves as a metal. On the other 
hand, in Figs. 82 and 83, a region where DOS has 
disappeared exists near the Fermi energy, and the 
system behaves as a Mott insulator. When densities of 
states corresponding to different values of L are 
remarked in Fig. 81, in case of L>4, the system has 
changed to an insulator. Therefore, in case of U=8, 
Mott metal-insulator transition has been confirmed to 
occur when L changes. 

In case of U^IO, the system behaves as a 
Mott insulator under all analyzed values of L, but 
since the value corresponding to a double of the value 
of 0 ) rendering DOS be zero is the gap as the Mott 
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insulator (Hubbard gap), it is appreciated that the 
width of the Hubbard gap increases with L. Therefore, 
it has been confirmed that Mott -Hubbard transition can 
be controlled by adjusting L. 

Figs. 84 through 86 show DOS of multiply- 
looped ring structures of a=0.5 modified by parallel 
movement proposed here. Setting U=8 in Fig. 84, U=10, 
in Fig. 85 and U=12 in Fig. numerical calculation was 
executed, selecting L=3, 4, 6, 8, 10, 20 for each value 
of U. This model is characterized in the tail 
structure of the density of state is smaller than the 
case of a=0 and may even disappears. That is, when a>0, 
the Hubbard gap width increases, and the tendency of 
the system as the Mott insulator is enhanced. This is 
the same as effect of parallel movement in a multiply- 
twisted helix (third embodiment), and it has been 
demonstrated that Mott metal -insulator transition in 
multiply- looped ring structures behave in the same 
manner as that in multiply- twisted helix. 

Next explained is a multiply- looped ring 
structure according to the 13th embodiment of the 
invention. In the 13th embodiment, quantum level 
statistics will be analyzed concerning quantum states 
in the multiply- looped ring structure to which random 
potentials are introduced, thereby to demonstrate that 
quantum chaos can be controlled by addition of magnetic 
impurities. 
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Explained below is an electron system on a 
multiply- twisted ring in the multiply-looped ring 
structure according to the 13th embodiment. 

Using a ring having L sites in one cycle, 
multiply- twisted rings are defined hierarchically. 
Assuming q hierarchies or layers, the total number of 
sites is Lq. The sites are described by a q- 
dimensional cube. 

a = (a^,a^,a^■aJ (228) 

where aic is an integer that satisfies I^Ok^L. 

Let the operator for generating a quantum at 
the a-th lattice point be Cj . Of course, let the 
operator satisfy the anticommutation relation 

(229) 

This quantum is a fermion having no free spin freedom. 
This corresponds to analysis of an electron in a solid 
in which spin orbit interaction can be disregarded. 

Here is used 

6(a-b)-Yls(ai-bj) (230) 

where 6(a-b) is a delta function defined by 

c/ fl when a=b 

b{a-b)=\ (231) 

10 otherwise 

Letting electrons be movable only among 
neighboring sites, determines their bonding 
relations. First introduced is the following function. 
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'E^^-\a,b;k,L) = 


(232) 


1 when b^-a^ = ±1 
. 1 when b^ -a,^ =±(1-L) 
0 otherwise 
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This function has been introduced to provide a 
periodical boundary condition and bonds between 
adjacent loops in the k-th layer. Next introduced is 


(233) 

j-i 

This is a function that takes 1 only when having the 
same suffix concerning layers beyond the k-th layer. 
Furthermore, by using 


Q{a,b;k,q) = 



- - i) 

1 


when 

when 


q >k +1 
q ^k +l 


(234) 


i?<*>(A:;a,L) = a^mod(r-l±— ,L)j (235) 

also introduced is 

Using these functions. 


is defined. When placing a=0, the above definition 
results in the same as that already discussed in the 
tenth embodiment. When a>0, the same effects as 
parallel movements of inter- layer bonds in the 
multiply- twisted helix is expected. Hereunder, using 
a=0.5, calculation will be conducted." 

Generally, a natural number j is defined as 
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(238) 


In this manner, sites of the multiply- looped ring 
structure can be equated with a one-dimensional lattice 

of l^j<Lq. The Hamiltonian H of this quantum system is 
5 defined as 

+ <239) 

a 

The second term of the Hcuniltonian is the term of the 
random potential. For each site, the random variable 
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15 


20 


--<v, <- (240) 

is generated to foim the Hamiltonian. The variable 
breadth v of the random potential is useful as a 
parameter determining the degree of the randomness. 
ta,b has been introduced for the random magnetic field. 
Here is determined as 

(241) 

where 0p,q=-0b,a is a random real niomber that satisfies 

(242) 

Due to hopping to the nearest -neighbor sites, the 

random phase factor 0a,b is re -suffixed, depending the 
sites. If the phase factor is integrated in the loop 
making one turn around a lattice point, a magnetic flux 
passing through the loop is obtained. This means a 
magnetic field is locally introduced to the random 
distribution of O<0a,b<23i this magnetic field is 
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absolutely random in both intensity and direction, and 
if it is averaged spatially, it becomes a zero magnetic 
field . 

When eigen energy of the Hamiltonian ^ is e , 
and the eigen vector is |m) , then 

= (243) 

where m=l, 2 , n. 

First, n quantum levels are standardized 
such that their nearest -neighbor level spacing becomes 
1 in average. The density of states of the system is 
defined by 

p(e)=-y<5(e-e„) (244) 

^ m»l 

and its staircase function 

^{^) = f_yvp{v) (245) 

is calculated. The staircase function obtained is 
converted by using a procedure called unfolding such 
that the density of states becomes constant in average. 
Using the quantum level obtained in this manner, the 
nearest -neighbor level spacing distribution P(s) and 

the A 3 statistics of Dyson and Mehta are calculated as 
quantum level statistics . As taught by Literatures 
(32) and (33), by using these statistics, it can be 
detected whether quantum chaos has been generated or 
not. It is also known that a quantum chaotic system is 
sensitive to perturbation from outside similarly to the 
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classical chaotic system, and analysis of quantum chaos 
is important as a polestar of designs of non-linear 
materials . 

In case of an integrable system, nearest - 
neighbor level spacing distribution P(s) and A 3 
statistics are those of Poisson Distribution 

Pp(5) = e‘* (246) 

^ 3 («) = ^ (247) 

In case of a quantum chaos under a magnetic field, it 
becomes GUE distribution 





(248) 


A3W 


1 


log(2n7i) + y 



+ 0(rt-‘) 


(249) 


where r is the Euler's constant. 

In this section, L, which determines the turn 
pitch of the multiply-looped ring structure, s, which 
determines the inter-layer bonding force, and v, which 
determines the intensity of the random potential, are 
varied as parameters. VIhen L= 6 , 10, 20 is used as the 
turn pitch, the values 


L 

q 

n 

6 

5 

7776 

10 

4 

10000 

20 

3 

8000- 
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were used as the total number of sites. Statistics 
values are calculated by using 2001 quantum levels, in 
total, taking 1000 states on each side from the band 
center among n states used for each L. 

Fig. 87 and Fig. 88 show quantum level 
statistics of multiply- looped ring structures of L =6 
while fixing the inter- layer bonding to s=l and 
changing the random potential intensity v. Fig. 89 and 
Fig. 90 show quanttim level statistics of multiply- 
looped ring structures of L=10 while fixing the inter- 
layer bonding to s=l and changing the random potential 
intensity v. Fig. 91 and Fig. 92 show quantum level 
statistics of multiply- looped ring structures of L=20 
while fixing the inter- layer bonding to s=l and 
changing the random potential intensity v. Figs. 87, 

89 and 91 show P(s) whilst Fig. 88 , 90 and 92 show A 3 
statistics values. In each of Figs. 87 through 92, 
when V is small, the system is in a metallic state, and 
the quantum level statistics is that of a quantum 
chaotic system of GUE. As v increases, the quantum 
level statistics changes toward the Poisson 
distribution. This is typical Anderson transition 
((34) Phys. Rev. 109, 1492 (1958); (35) Phys . Rev. Lett. 
42, 673(1979); (36) Rev. Mod. Phys. 57. 287(1985)). 

Through a careful review of the case of v=2, 
it is confirmed that the chaotic tendency is enhanced 
as L increases. On the other hand, in case of v=22, as 
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L increases, the system asymptotically changes toward 
the Poisson distribution. That is, as L increases, the 
level statistics values of this quantum system become 
more sensitive to changes of v. 

Similarly to the multiply- twisted helix, here 
is explained is a result of analysis of a system in 
which the inter- layer bonding intensity s is decreased. 
While fixing v=2. Fig. 93 shows A 3 statistics values of 
multiply-looped ring structures of L= 6 , Fig. 94 shows 
A 3 statistics values of multiply- looped ring structures 
of L=10, and Fig. 95 shows A 3 statistics values of 
multiply-looped ring structures of L=20. When s 
decreases from 1 , the inter-layer bonding becomes 
weaker, and the spatial dimension reduces. Therefore, 
it results in enhancing localization. It is 
appreciated from Figs. 93 through 95 that the system 
changes from a quantum chaotic sate to a considerably 
localized state along with changes from s=l to s= 0 . 4 . 
The degree of the changes becomes larger as L becomes 
smaller. This is because, when L is small, because of 
a small diameter of the basic double ring, this is 
substantially a one -dimensional tube, but in contrast, 
when L is large, because of a large diameter of the 
basic double ring, this itself forms a two-dimensional 
sheet . 

Next explained is a multiply- looped ring 
structure sys according to the 14th embodiment of the 
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invention. In the 14th embodiment, it is demonstrated 
that the nature of ferromagnetic transition occurring 
in a multiply- looped ring structure can be controlled 
by the turn pitch or parallel movement of inter- layer 
bonds . 

Explained below is a spin system on a 
multiply-twisted ring in the multiply-looped ring 
structure according to the 14th embodiment. 

Using a ring having L sites in one cycle, 
multiply- twisted rings are defined hierarchically. 
Assuming q hierarchies or layers, the total number of 
sites is Lq. The sites are described by a q- 
dimensional cube . 

a =(a„fl2>«3 •",) (250) 

where ak is an integer that satisfies l<ak^L. 

Hamiltonian H of this spin system is defined 
as follows. 

(251) 

Here is employed three-dimensional vectors 
having values at vertices of a regular octahedron. 

That is, Sa=(l,0,0), (-1, 0, 0), (0, 1, 0), (0, -1, 0), 

(0, 0, 1), (0, 0, -1) are allowable as the spin value. 
Letting spins be movable only among neighboring sites, 
defines their bonding relation. First introduced 
is the following function. 
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fl when b^-a^=±l 

1 when b^-a^ =±{1-L) 

0 otherwise 


This function has been introduced to provide a 
periodical boundary condition when loops of the period 
L are formed in the k-th layer. Next, 


k-i 


W(a,b;k)-n6(aj-bj) 

y-i 


(253) 


is a function that takes 1 only when having the same 

suffix concerning layers beyond the k-th layer. 8(a-b) 
is a delta function defined by 


6(a-b)=^\ 


1 when a = b 
0 otherwise 


(254) 


Furthermore, by using 

- j)6(b^^j - j) when q>k+l 


Q(a,b;k,q) = 


1 


R^~\x;a,L) = 6^ 


= 6 mod(r -1± — ,L) 

2 I 


when q +1 


aL 


(255) 


(256) 


also introduced is 

a is the quantity corresponding to parallel movement of 
inter-layer bonds in multiply- twisted helixes. What is 
analyzed here are multiply- looped ring structures of 
a=0 and multiply- looped ring structures of a=0.5 
modified by parallel movements. Using these functions. 


\,b ~ \,a 




(258) 
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is defined. Hereunder, s=l is assumed. 

For simplicity, in the subsequent procedures, 
a natural number j is defined as 

+ (259) 

In this manner, sites of the multiply-looped ring 
structure can be equated with a one -dimensional lattice 
of IsjsLq. 

For the purpose of calculating physical 
values at finite temperatures using statistical 
mechanics, a distribution function 

Z = (260) 

is introduced, where T is the temperature, and P=l/T 
and kB=l are used. Expected value concerning an 
arbitrary function /(sj of the spin variable is 
calculated by 

(261) 

When n is the total number of sites, expected value of 
spontaneous magnetization is 

M=-y(5.) (262) 

Spontaneous magnetization is calculated by 
using the Metropolis method that is one of the Monte 
Carlo methods. Using 100000 as the Monte Carlo steps, 
different turn pitches, L=3, 4, 6, 8, 10, 20, were 
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analyzed. For each of them, the following total 
numbers of sites were used. 


L 

q 

n 

3 

8 

6561 

4 

6 

4096 

6 

5 

7776 

8 

4 

4096 

10 

4 

10000 

20 

3 

8000 


Fig. 96 shows spontaneous magnetization in 
case of a=0 . Similarly to the ferromagnetic transition 
in multiply- twisted helixes, critical temperature 
apparently varies with L. Fig. 97 shows spontaneous 
magnetization in case of a=0.5. A shoulder-like 
structure observed in Fig. 96 of a large L disappears 
here, and changes of the phase transition temperature 
is remarkable . That is , excellent control has been 
realized. It has been thus confirmed that 
ferromagnetic transition in multiply- looped ring 
structures can be controlled by the turn pitch and 
parallel movements of inter-layer bonds. 

Next explained is a multiply-twisted helix 
according to the 15th embodiment of the invention. In 
the 15th embodiment, explanation will be made about a 
helix-based hierarchy of multiple connections in which 
bonds of s’' in multiply- twisted helixes are bonded with 
a force s and make up s’'*^ chains. As will be explained 
later, it will be demonstrated that, while this system 
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exhibits small dependency of ferromagnetic transition 
upon the number of turns N and temperature changes of 
the opposite sign from that of multiply- twisted helixes, 
it strongly depends on the value of s. Making use of 
this dependency, physical phenomena appearing in the 
helix-based hierarchy of multiple connections can be 
controlled. 

First, the multiply- twisted helix (MTH) is 
reviewed once again. Assuming a one -dimensional 
lattice, numbers are assigned as q= ", -1, 0, 1, 

The spin at the q-th site is designated by Sq. When 
Hamiltonian H of the spin system discussed below is 
defined by 

(263) 

P.9 

positional relation between spins is determined by 


employing 

as 

•^p.q 





1 

when 

It 

1 



- 

s 

when 

1 

II 





when 

mod(p,JV) = 0 

and 

q^p + N^ + f(N) 

J =J = . 


when 


and 

q=p + N^ + f(N^) 


s* 

when 

mod(p,A^") = 2 

and 

q=p + N'' + f(N^) 


0 

otherwise 





and a multiply- twisted helix is introduced. mod(a, b) 
is the remainder as a result of division of a by b. 
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The average of nearest -neighbor sites in a q- 
layered structure is obtained as 

/II 1 \ 

f (AT)==2 + 2+2 — + -Vt- - + -4- =4 + 2 — -- (265) 

’ [N iV’-'j N-1 

and when 

z„(iV) = 4 + 2/(iV-l) (266) 

is obtained. 

The helix-based hierarchy of multiple 
connections (HBHMC) explained in the 15th embodiment is 
defined as shown below. Assuming a one -dimensional 
lattice, numbers are assigned as x= •, - 1 . 0 , 1 , ". 

The spin variable at the x-th lattice point is denoted 
by S(x) . 

Hamiltonian Hi of a spiral is defined as 

follows . 

H^=-^S{x)-Six + l)-s'2S(x)S(x + N) (267) 

X X 

For the purpose of introducing high-order bonds, sites 
are added cximulatively . Ha describing secondary bonds 
is 

'2\^(x)S2ix,l) + S2ix-H)S(x + N^ + f{N))\ (268) 

mod(x,Ar)=0 

Hi describing tertiary bonds is 

^3 =-^ 2;[‘^W''^3W) + 53(r,l)-53(r,2) 

mod(j:^^)=l 

+ S^(x,2)S(x + N^ + f(N^))] (269) 

H 4 describing quaternary bonds is 
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^4 = + -^4 (-^>2) • S,(x,3) 

mod(jr,N^)-2 


+ S,{x,3)-S(x + N^ + f(N^)) ] 


Generally for q^2 



+ + +/(A^‘''')) ] (271) 

S', however, is conducted only with x that satisfies 
mod(x, N^‘^)=q-2. Using it, Hamiltonian defined by 


(272) 

is discussed. This site-to-site relation is shown in 
Fig. 98. 

As to the helix-based hierarchy of multiple 
connections defined above, here is calculated the 
expected value of the number of nearest -neighbor sites. 
Assume here that the main chain spiral is made up of M 
sites. Regarding this main chain spiral, the expected 
number of nearest -neighbor sites is the same as that of 
the multiply- twisted spiral. Assuming q layers, the 
number of sites added other than the main chain spiral 
is found to be 

num + ^ + - - + (273) 

This nuniq(N) can be calculated by using 

m-I 
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num^{N) = + (9 - ( 275 ) 

is obtained, which results in 

num^{N)^\\mnum(N)-^M—^^^ (lie) 

■' {N -\) 

Considering that all of the added sites each have two 
nearest -neighbor sites, the expected number of nearest- 
neighbor sites of the entire system is found to be 

^ \ (277) 

M +num^(N) N -N +1 

Fig. 99 shows z„(N) and 57„(iV) . What should be 
noted is the dependency on N and that, while z^(N) was 
a decreasing function, w„(iV) is an increasing function. 
In addition, since the deviation from 4 contains N*^, 
quickly asymptotically approaches 4 and loses 
dependency on N. From this phenomenon, in 
ferromagnetic transition expected in the helix-based 
hierarchy of multiple connections, it is expected that 
N dependency of the critical temperature oppositely 
behaves from that of a multiply- twisted helix, and the 
variation width is presumed to be small. However, this 
is about a system in which all bonds are equal in force, 
that is, s=l, and the following simulations demonstrate 
that the critical temperature can be favorably 
modulated by changing s . 

The helix-based hierarchy of multiple 
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connections analyzed here is made up of M=5000 sites. 
This spin system is analyzed by Mote Carlo simulation. 
Here is employed a three-dimensional vector having 
values at vertices of a regular octahedron. That is, 
( 1 , 0 , 0 ), (- 1 , 0 , 0 ), ( 0 , 1 , 0 ), ( 0 , - 1 , 0 ), ( 0 , 0 , 1 ), 
(0, 0, -1) are allowable as the spin value. For the 
purpose of calculating physical values at finite 
temperatures using statistical mechanics, a 
distribution function 



is introduced, where T is the temperature, and P=l/T 
and kB=l are used. Expected value concerning an 
arbitrary function f{Sj) of the spin variable is 
calculated by 

{f(Sj))-^'^f(Sj)e-^ (279) 

When n is the total number of sites, expected value of 
spontaneous magnetization is 

< 280 ) 

Spontaneous magnetization is calculated by using the 
Metropolis method that is one of the Monte Carlo 
methods. 100000 was used as the Monte Carlo steps, and 
N=3, 4, 6, 8, 10, 20 were analyzed as the turn pitch. 

Fig. 100 shows magnetization curves for 
various N when fixing s=l. It is first recognized that 
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dependency of the transition temperature for 
disappearance of spontaneous magnetization upon N is 
small. Then, a result has been obtained that the 
transition temperature increases as N increases, as 
assumed previously. 

Figs. 101 through 106 show magnetization 
curves for various s. Fig. 101 is for N=3, Fig. 102 is 
for N=4, Fig. 103 is for N=6, Fig. 104 is for N=8, Fig. 
105 is for N=10, and Fig. 106 is for N=20. In any of 
them, it is appreciated that the transition temperature 
is controlled well by changes of s. It has been 
confirmed that ferromagnetic materials having wide 
transition temperatures can be realized by setting the 
value of s . 

Next explained is a helix-based hierarchy of 
multiple connections according to the 16th embodiment 
of the invention. In the 16th embodiment, Anderson 
transitions in helix-based hierarchies of multiple 
connections are analyzed by using quantum level 
Tstatistics to demonstrate that wide quantum states from 
metallic states exhibiting quantum chaos to localized 
states can be realized. It is demonstrated that 
electron states do not depend on the number of turns N 
so much but strongly depends on the value of inter- 
layer bonds s. Making use of this dependency, physical 
phenomena appearing in helix-based hierarchies of 
multiple connections can be controlled. 
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The helix-based hierarchies of multiple 
connections explained in the 16th embodiment are 
defined as shown below. Assuming a one -dimensional 
lattice, numbers are assigned as x= --, -1, 0, 1, 

The operator for generating a quantum at the x-th 
lattice point is denoted by c^(r) . Let this operator 
satisfy the anticommutation relation 

(281) 

This quantum is a fermion having no free spin freedom. 
This corresponds to analysis of an electron in a solid 
in which spin orbit interaction can be disregarded. 

The Hamiltonian H-^ of spirals is defined as 

= -/^c^(r)c(r + 1) - 5^c^(r)c(r + TV) + ^ v^c^(x)c(r) + H.C. ( 282 ) 

XXX 

For the purpose of introducing high-order bonds, sites 

are added cumulatively. H 2 describing secondary bonds 
is 

H 2 (^’1) + ^2 (^4)c(r + AT ^ + / (JV^)) J+ H.C. (283) 

mod(x,A^)=0 

■^3 describing tertiary bonds is 

^3 = -5 2 [c\x)^2(x,l) + 4^ (r,l)4(j^,2) 

mo<l(x,A^ ^ )=1 

+ 4 (x,2)c(x + N^ + f(N^)) ]+ H.C. (284) 

H^ describing quaternary bonds is 

H^=-s 2 [c^ (j:)4 + 4 (^.1)4 (x,2) + 4 (^.2)4 (j:,3) 

mod(x,Ar2)=2 

+ f\x,3)c{x + NU f(N^)) ]+ H.C. ( 285 ) 
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- -s2't‘'w,w)+(f +1) 

\*"i 


+ ^l(x,q-l)c(x+N‘' + f(N^'^)) ]+H.C. (286) 

S', however, is conducted only with x that satisfies 
mod(x, N^'^)=q-2. Using it, the Hamiltonian defined by 


H = (287) 

will be discussed. 

The helix-based hierarchy of multiple 
connections analyzed here is made up of a main chain of 
5000 sites . Since it also includes other sites 
introduced to form inter-layer bonds, the total number 
of sites is large. The total number of sites M for 
each N is as shown below. 


N 

M 

3 

8684 

4 

7168 

6 

6150 

8 

5769 

10 

5573 

20 

5231 


Further introduced is 

/W = -int^|^j (288) 

where int(x) is the maximum integer that does not 
exceed x. The second term of the Hamiltonian is the 
20 ... term of the random potential. For each site, the 
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random variable 



V 

2 


( 289 ) 


is generated to form the Hamiltonian. The variable 
breadth v of the random potential is useful as a 
5 parameter determining the degree of the randomness. 

Here, N, which determines the turn pitch of multiply- 
twisted spirals, s, which determines the inter- layer 
bonding force, and v, which determines the intensity of 
the random potential, are varied as parameters. 

10 When eigen energy of the Hamiltonian // is , 

and the eigen vector is |/n) , then 

H\m) = eJ\m) (290) 

where m=l, 2, M. First, M quantum levels are 
standardized such that their nearest -neighbor level 
15 spacing becomes 1 in average. The density of states of 

the system is defined by 

= (291) 

and its staircase function 


X{e)= ^dr)p{rf) (292) 

20 is calculated. The staircase function obtained is 

converted by using a procedure called "unfolding" such 
that the density of states becomes constant in average. 
Using the quantum level obtained in this manner, the 
nearest -neighbor level spacing distribution P(s) and 
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the A 3 statistics of Dyson and Mehta are calculated as 
quantum level statistics. For the calculation of the 
quantum level statistics, here are used 1501 levels in 
total including 750 levels above and 750 levels below 
the band center. As taught by Literatures (32) and 
( 33 ), by using these statistics, it can be detected 
whether quantum chaos has been generated or not. It is 
also known that a quantum chaotic system is sensitive 
to perturbation from outside similarly to the classical 
chaotic system, and analysis of quantum chaos is 
important as a poles tar of designs of non-linear 
materials . 

In case of an integrable system, nearest - 
neighbor level spacing distribution P(s) and A 3 
statistics are those of Poisson Distribution 

P,{s) = e-^ (293) 

As(n) = ^ (294) 

If the system currently reviewed exhibits quantum chaos, 
-it becomes GOE distribution 

= (295) 




Jt 


log(2jm) + r-^-j 

8 4 




( 296 ) 


where T is the Euler's constant. 

Fig. 107 and Fig. 108 show quantum level 
statistics of helix-based hierarchies of multiple 
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connections of N=10 while fixing the inter-layer 
bonding to s=l and changing the random potential 
intensity v. Fig. 107 shows P(s) whilst Fig. 108 shows 
A 3 statistics values. In both statistics values, the 
system is in a quantum chaotic state when v is small, 
but as V increases, it changes to a localized state 
described by the Poisson distribution, and shows the 
existence of Anderson transition. 

In Figs. 109 and 110, calculation is 
conducted by using various values of N while fixing s=l 
and v=2. Fig. 109 shows P(s) while Fig. 110 shows A 3 
statistics values. Since the weak random potential of 
v =2 is introduced, a quantum-chaotic metallic state is 
observed. It is appreciated that this nature does not 
depend on N and the quantum level statistics 
substantially obeys the GOE distribution for each value 
of N. Thus it must change toward localization if v is 
increased, so it is appreciated that Anderson 
transition always exists in the above system and the 
dimensionality of the system is near three. 

However, it is possible to decrease the 
helix-based hierarchy of multiple connections by 
decreasing s.- Fig. Ill shows A 3 statistics values with 
helix-based hierarchies of multiple connections of N=10 
when changing s as s=l, 0 . 8 , 0 . 6 , 0 . 4 , 0.2 while fixing 
v=2. Although the random potential is as weak as v=2, 
as s decreases, A 3 asymptotically approaches that of 
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the Poisson distribution. In this kind of systems, 
since any stronger randomness merely results in 
enhancing localization, there is no Anderson transition. 
This is an evidence of a low-dimensionality. 

A review has been made about quantum level 
statistics with helix-based hierarchies of multiple 
connections having various N, s and v, and it has been 
confirmed that wide quantum states from metallic states 
exhibiting quantum chaos to localized states can be 
realized. 

Next explained is a helix-based hierarchy of 
multiple connections according to the 17th embodiment 
of the invention. In the 17th embodiment, Mott 
transitions in helix-based hierarchies of multiple 
connections are analyzed by using the Green's function 
method to demonstrate that wide quantum states from 
Mott insulator states to metallic states can be 
realized. It is demonstrated that electron states do 
not depend on the number of turns N so much but 
strongly depends on the value of inter-layer bonds s. 
Making use of this dependency, physical phenomena 
appearing in helix-based hierarchies of multiple 
connections can be controlled. 

The helix-based hierarchies of multiple 
connections explained in the 17th embodiment are 
defined as shown below. Assuming a one-dimensional 
lattice, numbers are assigned as x= •, -1, 0, 1, 
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The operator for generating an electron of the spin a 
at the x-th lattice point is denoted by c^(r) . Let this 
operator satisfy the anticonunutation relation 

(297) 

5 The Hamiltonian of spirals is defined as 

H, = -t'^cl(x)c„(x + l)-s'2cl(x)c^(x + N) , 

x,a x,a 

+ U^n^(x)n^(x) + H.C. (298) 

X 

where 

na(x) = cl(x)c„(x) (299) 

10 is the operator of the number of spin a electrons , and 

U is the electron-to-electron interaction intensity. 

For the purpose of introducing high-order bonds, sites 
are added cumulatively. H2 describing secondary bonds 
is 

15 ^2 =-5 y 

n»d(5>)-(krTri 


20 


+ U \n^^ix,l)n^^{x,l) + H.C. 

mod(^,^)=0 


(300) 


where 

7/3 describing tertiary bonds is 


2 2 [cl (x)$^„ (r,l) + 4% (^.1)4,0 (^>2) 

+ 4 ’„(x, 2 ) 4 (j: + Ar’ + /(Ar>)) ] 


(301) 
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where 


+ ^ + ]+//.C. (302) 

mod(?JW)=0 


10 


15 


«3.o /) = Ho ))^3,o » 

i ^4 describing quaternary bonds is 


(303) 


H, = -s 


where 


mod(x,A^^)-2a«T;i 

+ ^|,(r,3)c,(r + iVV/(^^)) ] 

+ ^ y[«4,t(^’l)\i(^>l) + «4.i(^’2)n^.i(^>^ ] 

mod(x,^)=0 

+ ^ (r,3)n^ j (r,3) ]+ H.C. 

^4,0 J) = H,c (x, j)H.c (x, j) 

Generally for q^2 

+iA^,g-'Oca(x + N‘^ + f(N‘’-^)) ] 

+ (x,k)n^^ (x,k) + H.C. 


(304) 


(305) 


(306) 


S', however, is conducted only with x that satisfies 
mod(x, N'*'^)=q-2. Then, 

was used. Discussed below is 


20 


H = 




(308) 


9=1 


that is formed by using the Hamiltonians of respective 
layers defined above. 
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The helix-based hierarchy of multiple 
connections analyzed here is made up of a main chain of 
5000 sites. Since it also includes other sites 
introduced to form inter-layer bonds, the total number 
of sites is large. The total number of sites M for 
each N is as shown below. 


N 

M 

3 

8684 

4 

7168 

6 

6150 

8 

5769 

10 

5573 

20 

5231 


Further introduced is 

/(x) = -int^^j (309) 

where int(x) is the maximum integer that does not 
exceed x. Here, N, which determines the turn pitch of 
multiply-twisted spirals, s, which determines the 
inter-layer bonding force, and U, which determines the 
intensity of electron- to-electron interaction, are 
varied as parameters . 

Here is made a review about half -filled of 
one electron per site. The density of states of this 
system can be calculated by using the same method as, 
example, the first embodiment. The density of 
states thus obtained is shown 'below, and based on it, 
discussion will be developed. Fig. 112 shows densities 
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of states in helix-based hierarchies of multiple 
connections of N=10 by fixing the layer-to-layer 
bonding as s=l and changing electron-to-electron 
interaction intensity U. (o=0 is the Fermi energy. If 
the density of states remains the system is in a 
metallic phase, and if it disappears, the system is in 
a Mott insulator phase. In case of U=6, the system is 
in a metallic phase, and over U=8, the systems are Mott 
insulators. This is typical Mott transition by changes 
of intensity of electron-to-electron interaction. In 
general, the value corresponding to a double the value 
of to where the density of states disappears is the 
Hubbard gap, and it is meant that, as it gets larger, 
correlation between electrons is intensive. 

In Fig. 113, densities of states were 
calculated for various values of N while fixing s=l and 
U=6. It is appreciated from the figure that the 
dependency of the densities of states upon N is small, 
the systems for all values of N exhibit a metallic 
phase. This small dependency on N is the unique nature 
of helix-based hierarchies of multiple connections, 
which is observed also in magnetic transition and 
Anderson transition. 

However, models exhibiting remarkable 
modulation under control of the value of s also 
characterize the helix-based hierarchies of multiple 
connections. Fig. 114 plots densities of states for 
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various s when fixing N=6 and U=6. Fig. 115 plots 
densities of states for various s when fixing N=10 and 
U=6. Fig. 116 plots densities of states for various s 
when fixing N=20 and U=6. It has been confirmed that, 
as s varies, the density of states under o)=0 changes 
from a finite to zero and Mott transition has occurred. 

A review has been made about Half -filled 
electron systems in helix-based hierarchies of multiple 
connections with various N, s and U, and it has been 
confirmed that wide electron systems from Mott 
insulators to Metallic states can be realized. 

Next explained is a helix- based hierarchy of 
multiple connections according to the 18th embodiment 
of the invention. In the 18th embodiment, critical 
helical spiral structures that are multiple 
hierarchical structures in which the number of sites 
does not depend on the number of turns and has 
generality will be explained. It will be demonstrated 
that, when layer -to -layer bonding is s=l, critical 
temperature Tc of ferromagnetic transition does not 
depend on N. On the other hand, since Tc is remarkably 
modified by s, making use of this dependency, physical 
phenomena appearing in the critical spiral structures 
can be controlled. 

There is a structure not depending on N 
between the multiply-twisted helix explained in the 
15th embodiment and the helix-based hierarchy of 
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multiple connections. This structure , i.e. the 
critical spiral structure, is defined below. Here is 
defined the Hamiltonian Hi of a spiral as follows. 

= -^S(r)-5(r + l)-5^5'(r)-5(r + A^) (310) 

X X 

For the purpose of introducing high-order bonds, sites 
are added cumulatively. Ha describing secondary bonds 
is 

H2=-s y[S(x)xj2(x) + a2(x) S(x + N'' + f(N)) ] (311) 

mod(i>)=.0 

Ha describing tertiary bonds is 

H,=-s 2[5(x)<73(r) + a3(r)-5(r + A^' + /(A^^)) ] (312) 

mod(jt, ^ )-l 

H 4 describing quaternary bonds is 

H,=-s '2JS(x)x7^(x) + a^(x)-S(x + N^ + f(N^)) ] (313) 

Generally for q^2 

+ + ] (314) 

3’, however, is conducted only with x that satisfies 
rood(x, N^'^)=q-2. Using it , Hamiltonian defined by 

" = 2 ( 315 ) 

9 = 1 

is discussed.. This site-to-site relation is shown in 
Fig. 117. 

Assume here that the main chain spiral is 
made up of M sites. Regarding this main chain spiral, 
the expected number of nearest -neighbor sites is the 
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same as that of the multiply- twisted helix. Assuming q 
layers, the number of sites added other than the main 
chain spiral is found to be 


(iV) = M + • • • + 

This NUMq(N) can be calculated by using 


r.2r-.^[r“‘-r] 

/n=l r - 1 


NUM(N) = -^fl- jV'-"! 


(318) 


is obtained, which results in 


NUM„ = UmNUM(N) = 

9 -** ’ 1 


(319) 


Considering that all of the added sites each have two 
nearest -neighbor sites, the expected number of nearest - 
neighbor sites of the entire system is found to be 

~ ^ 2NUM.(N) ]■ 4 (320) 

Therefore, the number of nearest -neighbor sites does 
not depend on N, but it is the general number 4. 
Therefore, freedom can be integrated by involving 
spiral structures of lower layers, the system can be 
regarded as a structure of a critical region, which the 
involved group acts on. However, this discussion is 
limited to the case where all bonds have the same 
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intensity, that is, in case of s=l, and the following 
simulation will demonstrates that, by changing s, the 
critical temperature can be modulated favorably. 

The critical spiral structure analyzed here 
5 is made up of Mo=7000 sites. The other number Mi of 

sites introduced to make layer- to -layer bonds for each 
N is shown below. 


N 

M 

3 

3483 

4 

2315 

6 

1377 

8 

977 

10 

752 

20 

331 
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15 


20 


Spins of the added sites described by o are numbered as 
Ok (k=l, 2, •" Ml) hereunder. 

This spin system is analyzed by Monte Carlo 
simulation. As spin freedom, a three-dimensional 
vector having values at vertices of a regular 
octahedron is employed. That is, (1,0,0), (-1, 0, 0), 
(07 1, 0), (0, -1, 0), (0, 0, 1), (0, 0, -1) are 


allowable as the spin value . 

For the purpose of calculating physical 
values at finite temperatures using statistical 
mechanics, a distribution function 


Z = \e 


PH 


( 321 ) 


is introduced, where T is the temperature, and p=l/T 
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and kB=l are used. The expected value of spontaneous 
magnetization is 
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1 Mo 1 M, 


(322) 


Spontaneous magnetization is calculated by using the 
Metropolis method that is one of the Monte Carlo 
methods. 100000 was used as the Monte Carlo steps, and 
N=3, 4, 6, 8, 10, 20 were analyzed as the turn pitch. 

As the function f(x) determining the bonds. 


/(r) = -int^^j (323) 

was used, int(x) is the maximum integer that does not 
exceed x. 

Fig. 118 shows magnetization curves for 

various N when fixing s=l. It is first appreciated 

that transition temperature for spontaneous 

magnetization disappears is the general value that does 

not depend on N. Behaviors near the critical 

temperature and magnetization curves depend on N, and 

it is suggested that the critical index in the critical 

region, for example, will depend on N. 

Figs. 119 through 124 show magnetization 

* 

curves for various s. N=3 is used in Fig. 119, N=4 in 
Fig. 120, N=6 in Fig. 121, N=8 in Fig. 122, N=10 in Fig. 
121, and N=20 in Fig. 122. In any of these cases, it 
is appreciated that transition temperature is favorably 
controlled by changes of s. It has been confirmed that 
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ferromagnetic materials having wide transition 
temperature can be realized by setting the value of s. 

Next explained is a critical spiral structure 
according to the 19th embodiment of the invention. In 
the 19th embodiment, Anderson transitions in critical 
spiral structures are analyzed by using quantum level 
statistics to demonstrate that wide quantum states from 
metallic states exhibiting quantum chaos to localized 
states can be realized. It is demonstrated that 
electron states do not depend on the number of turns N 
so much but strongly depends on the value of inter- 
Isycr bonds s. Making use of this dependency, physical 
phenomena appearing in critical spiral structures can 
be controlled. 

A quantum system in a critical spiral 
structure is defined. Assuming a one -dimensional 
lattice, numbers are assigned as x=*", -1, o, 1, •*•. 

The operator for generating a quantum at the x-th 
lattice point is denoted by c^(x) . Let this operator 
satisfy the anticommutation relation 

jc(r),c^(y) ( 324 ) 

This quantum is a fermion having no free spin freedom. 
This corresponds to analysis of an electron in a solid 
in which spin orbit interaction can be disregarded. 

The Hamiltonian of spirals is defined as 

^ X r 



For the purpose of introducing high-order bonds, sites 
are added cumulatively. Generally for qs2, 

'[c\x)d^(x) + dl(x)c(x + N^ + ]+H.C. (326) 

3', however, is conducted only with x that satisfies 
5 mod(x, N**'^)=q-2. The operator c^(r) used here is an 

operator for generating a quantum at a site added, and 
it satisfies 

( 327 ) 

and 

10 ^»,rf;(y)}=0 (328) 

Using it , the Hamiltonian defined by 

(329) 

will be discussed. 

The critical spiral structure analyzed here 
15 is made up of Mo=7000 sites. The other number Mi of 

sites introduced to make layer -to- layer bonds for each 
N is shown below. 
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20 Further employed was 




( 330 ) 
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int(x) is the maximum integer that does not exceed x. 
The second term of the Hamiltonian is the term of the 
random potential. For each site, the random variable 
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is generated to form the Hamiltonian. The variable 
breadth v of the random potential is useful as a 
parameter determining the degree of the randomness . 
Here, N, which determines the turn pitch of multiply- 
twisted spirals, s, which determines the inter-layer 
bonding force, and v, which determines the intensity of 
the random potential, are varied as parameters. 

When eigen energy of the Hamiltonian H ±s e . 
and the eigen vector is |»i) , then 

H\m) = eJ[m) (332) 

where m=l, 2, n (n=Mo"l‘Mi). First, n quantum levels 

£„ are standardized such that their nearest -neighbor 
level spacing becomes 1 in average. The density of 
states of the system is defined by 

2 n 

p(e)=-Yd(s-sJ (333) 

and its staircase function 


(334) 

is calculated. The staircase function obtained is 
converted by using a procedure called "unfolding" such 
that the density of states becomes constant in average. 
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Using the quantum level obtained in this manner, the 
nearest -neighbor level spacing distribution P(s) and 
the A 3 statistics of Dyson and Mehta are calculated as 
quantum level statistics. For the calculation of the 
5 quantum level statistics, here are used 1501 levels in 

total including 750 levels above and 750 levels below 
the band center. As taught by Literatures (32) and 
(33), by using these statistics, it can be detected 
whether quantum chaos has been generated or not . It is 

10 also known that a quantum chaotic system is sensitive 

to perturbation from outside similarly to the classical 

chaotic system, and analysis of quantum chaos is 

important as a polestar of designs of non-linear 
materials . 

15 In case of an integrable system, nearest - 

neighbor level spacing distribution P(s) and A 3 
statistics are those of Poisson Distribution 

Pp{s) = e~’ ( 335 ) 

A 3 (/i) = ^ (336) 

20 If the system currently reviewed exhibits quantum chaos, 

it becomes GOE distribution 

(337) 

1 5 ' 

^zin) = — \og{7joi) + Y - — +0{n~^) (338) 

Jl O 4 

where 7 is the Euler's constant. 
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Fig, 125 and Fig. 126 show quantum level 
statistics of critical spiral structures of N=10 while 
fixing the inter-layer bonding to s=l and changing the 
random potential intensity v. Fig. 125 shows P(s) 
whilst Fig. 126 shows A 3 statistics values. In both 
statistics values, the system is in a quantum chaotic 
state when v is small, but as v increases, it changes 
to a localized state described by the Poisson 
distribution, and shows the existence of Anderson 
transition. 

In Figs. 127 and 128, calculation is 
conducted by using various values of N while fixing s=l 

and v=2. Fig. 127 shows P(s) while Fig. 128 shows A 3 
statistics values. Since the weak random potential of 
v =2 is introduced, a quantiam- chaotic metallic state is 
observed. It is appreciated that this nature does not 
depend on N and the quantum level statistics 
substantially obeys the GOE distribution for each value 
of N. Thus it must change toward localization if v is 
Increased, so it is appreciated that Anderson 
transition always exists in the above system and the 
dimensionality of the system is near three. 

However, it is possible to decrease the 
critical spiral structure by decreasing s. Fig. 129 

shows A 3 statistics values with critical spiral 
structures of N=10 when changing s as s=l, 0.8, 0.6, 

0.4, 0.2 while fixing v=2. Although the random 



potential is as weak as v= 2 , as s decreases, A 3 
asymptotically approaches that of the Poisson 
distribution. In this kind of systems, since any 
stronger randomness merely results in enhancing 
localization, there is no Anderson transition. This is 
an evidence of a low-dimensionality. 

A review has been made about quantum level 
statistics with critical spiral structures having 
various N, s and v, and it has been confirmed that wide 
quantum states from metallic states exhibiting quantum 
chaos to localized states can be realized. 

Next explained is a critical spiral structure 
according to the 20th embodiment of the invention. In 
the 20th embodiment, Mott transitions in critical 
spiral structures are analyzed to demonstrate that wide 
quantum states from Mott insulator states to metallic 
states can be realized. It is demonstrated that 
olectron states do not depend on the number of turns N 
so much but strongly depends on the value of inter- 
layer bonds s. Making use of this dependency, physical 
phenomena appearing in critical spiral structures can 
be controlled. 

The electron system in a critical spiral 
structure is defined. Assuming a on -dimensional 
lattice, numbers are assigned as x= --, -1, 0, 1, •. 

The operator for generating an electron of the spin o 
at the x-th lattice point is denoted by c^(r). Let this 



( 339 ) 


operator satisfy the anticonunutation relation 

The Hamiltonian of spirals is defined as 

x,a x,o 

5 +U'^n^{x)n^{x)-\-H.C. (340) 

X 

where 

K(x)=cl{x)c^{x) (341) 

is the operator of the number of spin o electrons, and 
U is the electron-to-electron interaction intensity- 
10 For the purpose of introducing high-order 

bonds, sites are added cumulatively. Generally for qa2, 

+ Ul'm^^{x)m^^{x)-^H.C. (342) 

S', however, is conducted only with x that satisfies 
15 mod(x, N‘*‘^)=q-2. The operator cj<,(r) used here is an 

operator for generating an electron of spin a at a site 
added, and it satisfies 

(>')}= ^x,y^,,p^p ( 343 ) 

and 

20 ^p(^ldlo(y)\=0 (344) 

Additionally, the notation 

= (345) 

was used. Using them, the Hamiltonian defined by 

( 346 ) 

9-1 

25 will be discussed. 
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The critical spiral structure analyzed here 
is made up of a main chain with Mo=7000 sites. The 
other number Mi of sites introduced to make layer-to- 
layer bonds for each N is shown below. 
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Further employed was 

/(r) = -int|^j (347) 

int(x) is the maximum integer that does not exceed x. 
Here, N, which determines the turn pitch of multiply- 
twisted spirals, s, which determines the inter-layer 
bonding force, and v, which determines the intensity of 
the random potential, are varied as parameters. 

Here is made a review about half-filled of 
one electron per site. The density of states of this 
system can be calculated by using the same method as, 
for example, the first embodiment. The density of 
states thus obtained is shown below, and based on it, 
discussion will be developed. Fig. 130 shows densities 
of states in critical spiral structures of N=10 by 
fixing the layer -to -layer bonding as s=l and changing 
electron-to-electron interaction intensity U. oj=0 is 



the Fermi energy. If the density of states remains the 
system is in a metallic phase, and if it disappears, 
the system is in a Mott insulator phase. In case of 
U=6, the system is in a metallic phase, and over U=8, 
the systems are Mott insulators. This is typical Mott 
transition by changes of intensity of electron-to- 
electron interaction. In general, the value 
corresponding to a double the value of co where the 
density of states disappears is the Hubbard gap, and it 
is meant that, as it gets larger, correlation between 
electrons is intensive. 

In Fig. 131, densities of states were 
calculated for various values of N while fixing s=l and 
U=6. It is appreciated from the figure that the 
dependency of the densities of states upon N is small, 
the systems for all values of N exhibit a metallic 
phase. This small dependency on N is the unique nature 
of critical spiral structures, which is observed also 
in magnetic transition and Anderson transition. 

However, models exhibiting remarkable 
modulation under control of the value of s also 
characterize the critical spiral structures. Fig. 132 
plots densities of states for various s when fixing N=6 
and U=6. Fig. 133 plots densities of states for 
various s when fixing N=10 and U=6. Fig. 134 plots 
densities of states for various s when fixing N=20 and 
U=6. It has been confirmed that, as s varies, the 



density of states under co=0 changes from a finite to 
zero and Mott transition has occurred. 

A review has been made about Half -filled 
electron systems in helix-based hierarchies of multiple 
connections with various N, s and U. and it has been 
confirmed that wide electron systems from Mott 
insulators to Metallic states can be realized. 

Next explained is the 21st embodiment of the 
invention. In the 21st embodiment, a manufacturing 
method of a multiply- looped ring structure of q=2 and 
N=8, taken as an example, will be explained. Assume 
that the multiply- looped ring structure is explained by 
using a cylindrical quantum box by compound 
semiconductor hetero- junction. Essentially, the same 
manufacturing method is applicable also to multiply- 
twisted helices . 

That is, in the 21st embodiment, as shown in 
Fig. 135, a GaAs layer 2 is epitaxially grown on an 
AlGaAs substrate 1 by, for example, metal-organic 
chemical vapor deposition (MOCVD) . The GaAs layer 2 is 
grown up to, for example, 10 nm. 

After that, a Si02 film 3 is formed on the 
GaAs layer 2 by, for example, vacuum evaporation or 
chemical vapor deposition (CVD) . Further formed 
thereon is a circular resist pattern (not shown) in a 
ring-shaped arrangement as a whole by, for example, 
electron beam lithography. Thereafter, using the 
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resist pattern as a mask, the Si02 film 3 is 
selectively removed by, for example, reactive ion 
etching (RIE) . The resist pattern is removed 
thereafter. In this manner, the circular Si02 film 3 
is formed on the GaAs layer 2 in a ring-shaped 
arrangement as a whole. 

Next as shown in Fig. 136, using the Si02 
film 3, the GaAs layer is etched in a direction noinnal 
to the substrate surface by, for example, RIE. In the 
etching process by RIE, a chlorine-based gas is used as 
the reaction gas (etching gas ) . This is for the 
purpose of using the AlGaAs substrate 1 as the etching 
stop layer to ensure that the GaAs layer 2 is 
selectively etched. 

Next as shown in Fig. 137, using the Si02 
film 3 as a growth mask, an AlGaAs layer 4 is 
epitaxially grown by, for example, MOCVD to bury 
portions of the GeiAs layer 2 removed by the etching 
with the AlGaAs layer 4 . 

Then, after removing the SiOa film 3 by 
etching using a hydrofluoric acid, an AlGaAs layer 5 is 
epitaxially grown on the entire surface as shown in Fig. 
138 by, for example, MOCVD. The AlGaAs layer 5 is 
grown up to the thickness of, for exeunple, 5 nm. 

Through these steps, a quantum box by AlGaAs/GaAs 
hetero- junction, which is the element forming the 
multiply-looped ring structure, is obtained. 
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Those steps may be repeated while rotating 
the pattern of the SiOa film 3 of Fig. 139 by a 
predetermined value every time to finally obtain the 
multiply-looped ring structure. 

Molecular beam epitaxy (MBE) is also 
applicable as the growth method other than MOCVD. The 
shape of the quantum box need not be cylindrical, but 
any other appropriate shape, such as a quadratic prism, 
may be used. 

Heretofore, embodiments of the invention have 
been explained specifically; however, the invention is 
not limited to those embodiments but can be modified in 
various modes based on the technical concept of the 
invention. 

As described above, according to the 
invention, it is possible to obtain new highly- 
functional materials that are complementary with 
fractal -shaped materials and exhibit new physical 
properties by making up a one-dimensional structure, 
spiral structure or ring structure from at least two 
layers of thinner one-dimension unit structures, spiral 
unit structure or ring unit structures bonded to each 
other in at least one site and setting the curvature, 
turn pitch or number of elements thereof such that the 
structure exhibits a nature regulated thereby, or by 
making up a one- dimensional structure, spiral structure 
or ring structure from thinner one -dimensional unit 
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structures, spiral unit structures or ring unit 
structures and setting the curvature, turn pitch or 
number of elements thereof such that the structure has 
a dimensionality regulated thereby. 
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